Historic,  Archive  Document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


SP^f  FERE  DISTANCE 
FROM  BURNING  TREES 

A  PREDICTIVE  MODEL 

FRANK  A.  ALBINI 


■V- 


V  > 


ij 


^SDA^f^tional  Agriculture,,;,^ 

10301  Baltimore  Blvd 


USDA  Forest  Service 
General  Technical  Report  INT-56 

July  1979 


SPOT  FIRE  DISTANCE  FROM  BURNING  TREES  — 

A  PREDICTIVE  MODEL 

F.  A.  Albini 


USDA  National  Agricultural  Library 
NAL  Building 
10301  Baltimore  Blvd. 
Beltsville,  MO  20705-2351 


INTERMOUNTAIN  FOREST  AND  RANGE  EXPERIMENT  STATION 

Forest  Service 
U.S.  Department  of  Agriculture 
Ogden,  Utah  84401 


THE  AUTHOR 


FRANK  A.  ALBINI  is  a  Mechanical  Engineer,  assigned  to  the  Fire 
Fundamentals  research  work  unit  at  the  Northern  Forest  Fire 
Laboratory  in  Misssoula,  Montana.     He  earned  a  Ph.D.  from  the 
California  Institute  of  Technology  in  1962,  where  he  also  ob- 
tained his  undergraduate  training   (B.S.  1958,  M.S.  1959).  He 
joined  the  Forest  Service  in  October  1973  after  12  years  of 
pure  and  applied  research  and  systems  analysis  both  in  pri- 
vate industry  and  at  the  nonprofit  Institute  for  Defense 
analyses . 


ACKNOWLEDGMENTS 


This  model  incorporates  a  cross  section  of  the  technical 
fields  associated  with  wildland  fire  science.     The  author  is  in- 
debted to  specialists  from  several  disciplines  for  their  assist- 
ance and  for  educated  guesswork  needed  to  fill  in  pieces  of  the 
models  presented  here.     The  author  accepts  full  responsibility 
for  any  misapplication  of  data  or  misinterpretation  of  theories, 
but  wishes  to  acknowledge  these  sources: 

Dr.  James  K.  Brown  -  for  data  on  foliage  weights  and  sur- 
face/volume ratios  and  introduction  to  the  literature  in  which 
such  data  reside. 

Dr.  Ralph  Wilson,  Jr.,   -  for  data  on  firebrand  burning  rate 
and  consultation  that  helped  in  understanding  the  data. 

Dr.  Donald  Fuquay,  Mr.  Robert  Baughman,  and  Dr.  Donald 
Latham  -  for  encouragement,  guidance,  and  constructive  criticism 
in  dealing  with  the  poorly-defined  meteorological  aspects  of  thi 
problem. 

All  of  these  people  are  located  at  the  Intermountain  Sta- 
tion's Northern  Forest  Fire  Laboratory  in  Missoula,  Montana. 


RESEARCH  SUMMARY 


This  paper  presents  a  predictive  model  for  calculating  the 
maximum  spot  fire  distance  to  be  expected  when  firebrands  are 
thrown  into  the  air  by  the  burning  of  tree  crowns  either  individ- 
ually or  in  small  groups.     Variables  included  in  the  model  are 
the  quantity  and  surf ace/volume  ratio  of  foliage  in  the  burning 
tree(s),  the  height  of  the  tree(s),  and  the  wind  field  that  trans- 
ports the  firebrands,  and  the  firebrand  burning  rate.     Many  aspects 
of  the  processes  modeled  are  as  yet  incompletely  understood  and  so 
are  only  weakly  represented.     Improved  models  are  needed  for  all 
of  the  processes  before  predictions  can  be  made  v/ith  certainty, 
but  a  step-by-step,  graphical  procedure  is  presented  here  for  test 
and  evaluation  in  the  field.     No  validation  data  are  available  at 
present. 
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INTRODUCTION 


In  tlie  jargon  of  forest  fire  control  a  '"spot  fire"  is  a  fire  set  outside  tlie 
perimeter  of  the  main  fire  by  flying  sparks  or  embers   (USDA  Forest  Service  1956} .  A 
wildland  fire  is  said  to  be  "spotting"  when  it  produces  sparks  or  embers  that  are 
carried  by  the  wind  aiid  start  nev;  fires  beyond  the  zone  of 'direct  ignition  by  the 
main  fire.     This  type  of  fire  spread  occurs  with  increasing  frequency  as  fuel  mois- 
ture content  declines   (table  1) . 

Fire  spread  by  spotting  is  a  chance  event.     'Vlien  it  will  occur  and  over  what 
distance  are  at  present  unpredictable  quantities  that  complicate  the  task  of  wild- 
fire control.     And  fires  ignited  for  desirable  effects  under  prescribed  conditions 
sometimes  escape  their  intended  boundaries  because  of  spotting. 

This  vexing  phenomenon  has  long  challenged  both  the  research  community  and 
fire  management  practitioners.     This  paper  presents  a  theoretical  framework  and 
instructional  materials  for  predicting  the  maximum  spot  fire  distance  from  burning 
trees.     With  testing  and  refinement,  the  model  given  here  should  enable  researchers 
and  managers  alike  to  better  understand  and  cope  with  spot  fires . 


Problem  Description 

Spotting  occurs  over  a  wide  range  of  distances,  depending  upon  windspeed  and 
the  type  (and  quantity)  of  fuel  involved  in  the  main  fire  as  well  as  fuel  moisture. 
Tlie  severity  of  a  potential  spotting  problem  can  be  described  numerically  by  the 
maximum  spot  fire  distance  to  be  anticipated  under  the  conditions  in  question.  The 
problem  addressed  in  this  paper  arises  under  conditions  of  intermediate  fire  severity 
in  which  spotting  distances  up  to  a  mile  or  two  might  be  encountered.     These  burning 
conditions  represent  a  rougli  upper  limit  to  fire  severity  under  which  prescribed 
burning  might  be  carried  out.     In  some  cases  an  estimate  of  maximum  possible  spot- 
ting distance  may  assist  in  delimiting  prescription  regimes.     Large  wildfires  also 
occur  frequently  under  burning  conditions  that  admit  intermediate -range  spotting. 
This  situation  arises  when  fire  intensity  is  not  so  high  as  to  cause  extensive 
crowning  or  crown  fire  spread  (Van  Wagner  1977} ,  but  when  the  intensity  is  suffi- 
ciently high  that  the  fire  cannot  be  "trusted"  to  burn  only  surface  fuels  (table  1} . 

Situations  not  oonsideved  in  this  work  are  those  extreme  cases  in  which  spotting 
may  occur  up  to  tens  of  miles  from  the  main  front,  as  in  running  crown  fires,  fires 
in  heavy  slash  or  chaparral  under  extreme  winds,  and  fires  in  which  fire  whirls  loft 
burning  material  high  into  the  air. 

The  situation  treated  here  starts  with  the  assumption  that  there  occurs  the 
occasional  or  sporadic  "torching  out"  of  individual  trees  or  small  groups  of  trees. 
These  singular  events  provide  the  mechanism  for  the  lofting  of  firebrands  that  are 
then  carried  by  the  prevailing  wind.     IVhile  this  is  not  the  most  severe  case  of 
spread  by  spotting,  it  is  a  situation  that  often  confronts  planners  in  fire  sup- 
pression and  fire  prescription  work. 

A  method  is  sought  that  will  permit  prediction  of  the  maximum  spot  fire  dis- 
tance to  be  anticipated  from  burning  trees,  given  a  knowledge  of  the  terrain,  forest 
cover,  and  windspeed.     If  maximum  spotting  distance  can  be  predicted,  the  information 
might  be  useful  in  prescription  writing,  selecting  fire  control  line  placement, 
positioning  of  "spot  chasers,"  and  similar  matters  that  fire  managers  must  deal  with. 
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Table  I .  --Spotting  potential  related  to  fuel  moisture  content 

(This  table  is  based  on  data  compiled  by  Hal  E.  Anderson,  USDA  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station,  Northern  Forest  Fire  Labora- 
tory, for  Fire  Behavior  Officer  training  ('National  Interagency  Fire  Training 
Center,  Marana  Air  Park,  Arizona).)  s 


Fuel  moisture  content  Relative  ease  of  chance  ignition,  likelihood 

Relative  1/2-inch  Forest        of  spotting,  general  burning  conditions 

humidity*  fuel  stick  litter 
______    Percent  -_---- 


>25 

Little  or  no  ignition-' 

>60 

>15 

>20 

Very  little  ignition;^  some  spotting  may  occur 
with  winds  above  9  mi/h^ 

45-60 

12-15 

15-19 

Low  ignition  hazard  -  campfires  become 
dangerous;^  glowing  brands  cause  ignition 
when  relative  humidity  <50  percent^ 

40+ 
30-45 

7-12 

11-14 

Medium  ignitibility  -  matches  become 
dangerous;-'  "easy"  burning  conditions'* 

26-40 

8-10 

High  ignition  hazard  -  matches  always 
dangerous;-'  occasional  crowning,  spotting 
caused  by  gusty  winds;^  "moderate"  burning 
conditions^ 

15-30 

5-7 

5-7 

Quick  ignition,  rapid  buildup,  extensive 
crowning;  any  increase  in  wind  causes  increased 
spotting,  crowning,  loss  of  control;^  fire 
moves  up  bark  of  trees  igniting  aerial  fuels; 
long  distance  spotting  in  pine  stands;^ 
dangerous  burning  conditions'* 

<15 

<5 

<5 

All  sources  of  ignition  dangerous;'  aggressive 
burning,  spot  fires  occur  often  and  spread 
rapidly,  extreme  fire  behavior  probable;^ 
critical  burning  conditions'* 

*Relative 

humidity  is 

a  surrogate 

for  "fine"  fuel  moisture  content. 

^Gisbome  (1936) 

^USDA  Forest  Serv. ,  Northern  Region  (1973) 
^Florida  Div.  For.  (1973) 
'*Barrows  (1951) 


Approach  to  Solution 

The  approach  used  in  constructing  this  model  is  to  examine  each  phase  of  the 
process  and  develop  a  mathematical  description  for  it.     The  individual  "submodels" 
are  based  on  physical  principles  to  a  large  extent,  but  assumptions,  approximations, 
and  inadequately  supported  empirical  relationships  are  sprinkled  throughout.  Each 
submodel  provides  an  element  of  information  needed  to  predict  the  maximum  spot  fire 
distance.     As  better  theories  and  data  become  available,  the  submodels  can  be  re- 
vised and  improved.     But  the  framework  that  links  each  of  the  submodels  rests  on  a 
series  of  basic  assumptions  that  form  a  conceptual  picture  of  the  overall  process. 
The  process  modeled  is  the  following  sequence; 

A  tree,  or  a  small  group  of  trees,  "torches  out."    The  flame,  and  the  buoyant 
plume  above  it,  exist  for  a  brief  period  of  time.     This  fluid  flow  field  is  capable 
of  lofting  potential  firebrands  into  the  air.     The  flow  structures  are  described  by 
separate,  steady-state  models,  joined  at  the  tip  of  the  flame.     These  models  assume 
still  ambient  air  for  simplicity,  and  to  ensure  that  a  maximum  height  is  predicted 
for  the  firebrand  particles. 
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The  flight  of  an  inert  wood  c/linder  in  the  flame/plume  flow  field  serves  to 
predict  the  height  of  a  potential  firebrand  as  a  function  of  time.     The  particle  is 
assumed  to  start  from  the  top  of  a  tree  involved  in  the  flame  at  the  instant  the 
flame  is  established.     It  travels  vertically  until  the  flame  goes  out  and  the  flow 
field  collapses.      The  height  achieved  by  the  particle  depends  on  its  diameter  and 
density  and  on  the  flow  field  description. 

When  the  particle  reaches  its  maximum  height,  it  is  then  subjected  to  the  am- 
bient wind  field  and  transported  laterally  as  it  falls.     The  burning  of  the  particle 
is  accounted  for  by  a  simple  empirical  model  that  describes  its  falling  velocity  as 
a  function  of  time.     The  maximum  spot  fire  distance  is  established  by  requiring  that 
the  particle  is  totally  consumed  just  as  it  returns  to  the  ground.     Smaller  particles 
would  travel  farther  but  burn  out  before  reaching  the  ground;  larger  ones  could  not 
travel  so  far. 

The  wind  field  that  carries  the  firebrand  can  be  modeled  by  a  logarithmic  varia- 
tion with  height  for  flat  terrain  (fig.  1).     But  another  submodel  is  required  to 
reflect  the  influence  of  ridges  and  valleys  on  the  wind  field  (fig,  2). 

Each  of  the  submodels  is  discussed  more  fully  in  the  section  entitled  Synopsis 
of  Submodel  Content  (p.  19)  and  is  described  in  appendixes  A-F. 
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Factors  Omitted  in  this  Approach 


Before  this  model  can  be  applied  and  tested,  it  is  essential  that  the  user 
recognize  what  the  model  does  not  include.     Several  essential  features  of  the  spot- 
ting process  are  omitted  from  consideration  here  and  no  amount  of  adjustment  of  the 
submodels  will  allow  them  to  be  reflected  in  the  results-      These  are  factors  that 
must  be  incorporated  by  further  model  development.     Until  this  is  accomplished,  the 
user  must  supply  the  missing  pieces  by  whatever  means  he  can.     The  principal  missing 
elements  are: 


1.     The  likelihood  of  trees  burning.     No  consideration  is  included  here  of  the 
probability  that  a  tree  (or  group  of  trees)  will  "torch  out."    The  type  and  quantity 
of  surface  fuels,  the  burning  conditions,  overstory  species,  crown  separation,  etc., 
all  influence  this  factor.     The  model  presented  here  must  be  thought  of  as  a  "what 
if"  computation  aid;  tlie  user  must  supply  the  description  of  tlie  trees  that  are  to  be 
considered  to  burn  as  input  to  the  model.     Van  Wagner  (1977)  describes  conditions  for 
the  start  and  spread  of  crown  fires.     IVhen  burning  conditions  are  so  severe  that 
crown  fires  are  to  be  expected,  this  model  would  be  of  only  marginal  interest.  So 
Van  Wagner's  work  can  be  used  as  a  guide  to  conditions  under  which  one  should  consi- 
der the  chance  of  torcliing  and  spotting  to  be  significant. 
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2.  Availahility  of  optimum  firebrand  material.     This  model  presumes  that  at 
least  one  ideally  suited  firebrand  particle  exists  in  an  unfettered  state  near  the 
top  of  a  burning  tree.     This  is  consistent  with  the  intention  to  estimate  the  maxi- 
mum potential  spot  fire  distance  only.     Any  reduction  in  maximum  spotting  distance 
due  to  the  fact  that  an  ideal  particle  was  not  present  is  outside  the  scope  of  the 
model  presented  here.     The  fact  that  a  wood  cylinder  is  used  to  represent  the  maxi- 
mum-range firebrand  is  not  felt  to  be  a  source  of  significant  distortion.     This  as- 
sumption is  discussed  more  fully  below. 

3.  The  prob ability  of  spot  fire  ignition.     In  this  model  the  maximum-range 
particle  is  idealized  as  just  being  totally  consumed  at  the  instant  it  reaches  the 
ground.     The  fact  that  an  extremely  small  spark  or  ember  can  indeed  ignite  forest 
fuels  is  established  by  experimental  data,^  so  this  idealization  does  not  substan- 
tially weaken  the  model.     But  a  firebrand  must  come  into  contact  with  easily  ig- 
nited dry  fuel  (e.g.^  litter,  duff,  rotten  wood)  for  a  spot  fire  to  start.  This 
model  does  not  deal  with  the  chance  of  such  contact  or  the  probability  that  ignition 
will  occur  if  the  contact  is  made.     Here  we  present  only  the  distance  that  a  fire- 
brand can  travel  and  still  retain  the  possibility  of  starting  a  fire. 

4.  The  number  of  spot  fires.     The  probabilistic  question  of  how  many  spot 
fires  there  might  be  under  certain  conditions,  or  the  number  per  unit  surface  area, 
is  not  addressed  here.     Only  when  the  three  elements  listed  above  are  appropriately 
modeled  will  it  be  possible  to  predict  the  density  distribution  of  spot  fire  igni- 
tions  Qluraszew  and  Fedele  1976), 


ESTIMATING  MAXIMUM  SPOT  FIRE  DISTANCE 


This  section  gives  a  step-by-step  procedure  for  estimating  maximum  spot  fire 
distance  from  a  tree  {ot  group  of  trees)  that  "torches  out."    The  steps  are  numbered 
to  correspond  with  the  specimen  worksheet  (exhibit  1,  p.  8). 

1.     First,  record  the  species  of  tree  assumed  to  bum,  the  diameter  at  breast 
height  Cd.b.h.)  and  height  of  the  tree,  and  the  number  of  trees  assumed  to  bum 
simultaneously. 

These  data  will  be  needed  during  the  completion  of  the  worksheet,  and  it  is  im- 
portant that  the  investigator  "leave  some  tracks"  so  the  work  can  be  checked  or  con- 
tinued after  interruption.     Be  sure  to  cross  out  the  incorrect  or  circle  the  correct 
units  of  measurement  when  recording  dimensional  data.     With  few  exceptions  one  is 
best  advised  to  retain  one  system  of  units  for  each  worksheet  used. 


R.  A.  Wilson,  Jr.,  and  A.  P.  Brackebusch,  personal  coTranunicat ion  of  unpublished 
data.;  Northern  Forest  Fire  Laboratory,  Missoula,  Mont.,  September  1977. 
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2.  Consult  the  guide  in  table  2  and  select  from  figures  3-5  the  curves  to  use. 
Read  from  the  appropriate  curves  the  steady  flame  height  and  steady  flame  duration 
for  a  single  burning  tree  of  the  size  and  species  of  interest.     Record  these  on  the 
worksheet.     These  numbers  must  be  revised  if  more  than  a  single  tree  is  involved  in 
producing  the  flame. 

3.  Consult  figure  6  and  record  the  multiplying  factors  for  flame  height  and 
duration.     Note  that  the  flanie  height  multiplier  is  greater  than  unity  and  the  dura- 
tion factor  is  less  than  unity,  so  these  numbers  should  not  be  confused.     Carry  out 
the  indicated  multiplications  and  record  the  results.     It  is  not  necessary  to  record 
more  than  two  significant  figures.     For  ejcample,  a  duration  calculated  to  be  3.622 
should  be  recoided  as  3.6;  a  flame  height  of  129.4  should  be  recorded  as  130. 

4.  Divide  the  height  of  the  burning  tree  (or  the  average  height  of  the  burning 
group  of  trees)  by  the  adjusted  flame  height  determined  in  step  3.     Be  sure  to  use 
the  same  units  of  measurement  so  the  ratio  will  be  dimensionless .     Record  the  result 
and  use  it  to  select  the  proper  curve  in  figure  7.     The  dimensionless  flame  duration 
derived  in  step  3  is  used  to  enter  figure  7;  read  the  appropriate  curve  and  record 
the  lofted  firebrand  height/steady  flame  height  value.     Multiply  this  quantity  by 
the  adjusted  steady  flame  height  from  step  3  and  record  the  result.     To  this  result 
add  one-half  the  tree  height  to  obtain  the  initial  firebrand  height  above  ground. 
This  result  is  used  in  figure  8  (step  6). 

5.  The  graphs  for  calculating  spotting  distance  use  the  windspeed  at  treetop 
height  at  the  average  elevation  (e.g.,  midslope  elevation  for  a  ridge-valley  system). 
It  is  often  necessary  to  adjust  the  available  data  on  windspeed  to  infer  this  value. 
First,  record  for  reference  the  elevation  and  height  above  terrain  at  the  spot  the 
windspeed  value  is  known  or  estimated.     Record  also  this  reference  windspeed,  the 
average  elevation,  and  the  mean  treetop  height  in  the  area  downwind  of  the  firebrand 
source.     This  mean  treetop  height  is  intended  to  characterize  the  general  forest 
cover  of  the  terrain  as  it  influences  the  wind  field  that  will  transport  a  fire- 
brand.    If  the  area  has  broken  forest  cover,  use  half  the  treetop  height  of  the 
forest-covered  portion. 

Subtract  the  average  terrain  elevation  from  the  elevation  where  the  wind  was 
measured.     Use  this  difference  in  figure  2  to  estimate  the  ratio  of  windspeed  at 
the  observation  point  to  that  at  the  average  elevation.     Record  this  value  in  the 
blank  space  marked  X  on  the  worksheet. 

Next,  a  correction  may  be  needed  for  the  height  above  the  vegetation  at  which 
the  reference  windspeed  was  determined.     The  reference  windspeed  should  be  reduced 
if  it  was  measured  above  the  canopy  top  represented  by  the  mean  treetop  height  noted 
above.     If  so,  divide  the  windspeed  measurement  height  by  the  mean  treetop  height 
and  obtain  from  figure  1  the  windspeed  ratio.     Record  this  value  in  the  blank  space 
marked  Y. 


^Figures  are  given  for  three  regions:     Intermountain  west,  north  central/north- 
east, and  south/southeast.     The  literature  revealed  no  foliage  weights  for  west/ 
northwest  tree  species.     Intermountain  species  might  serve  as  adequate  substitutes, 
such  as  grand  fir  for  noble,  red,  and  Pacific  silver  firs,  subalpine  fir  for  white 
fir,  western  white  pine  for  sugar  pine  and  Monterey  pine,  and  Ponderosa  pine  for 
Jeffrey,  Coulter,  and  Digger  pine.     Species  that  are  essentially  the  same,  such  as 
Ponderosa  pine,  Douglas-fir,  or  lodgepole  pine,  should  be  directly  substituted  until 
more  data  become  available. 
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If  the  u'indspeed  reference  height  is  not  greater  than  the  mean  treetop  height, 
then  this  correction  is  not  to  be  made  using  figure  1.     Because  no  generally  ac- 
cepted procedure  is  known  to  the  author  for  correcting  windspeed  measurements  icp-jard 
to  treetop  height ,  the  following  suggestions  are  offered  only  as  interim  measures 
and  are  not  based  upon  accepted  theory  or  specific  data.^ 


Windspeed  Measurement  Site  Y  Valve  Suggested 

Open  ridgetop,  20-foot  anemometer  tower  1.0 
Open  ridgetop,  hand-held  anemometer  .5 
Forest  opening,  20-foot  anemometer  tower  .5 
Forest  opening,  hand-held  anemometer  .4 
Under  canopy,  hand-held  anemometer  .3 


Apply  the  corrections  for  elevation  (X)  and  height  above  terrain  (YJ  by  using 
these  factors  as  divisors  of  the  reference  windspeed  value.     Enter  the  adjusted 
windspeed  value  on  the  worksheet.     This  is  used  in  the  next  step. 

6.  Figure  8  is  a  nomograph  allowing  one  to  predict  maximum  spot  fire  distance 
over  flat  terrain  with  a  uniform  forest  cover.     This  distance  can  be  corrected  for 
the  effects  of  terrain  on  the  wind  field,  so  must  be  computed  even  if  the  terrain  is 
not  approximately  flat.     T\vo  versions  of  figure  8  are  given- -8A  for  British  units 

of  measurement  and  8B  for  metric  units.     Qioose  the  appropriate  version  of  the 
figure  and  proceed  as  follows: 

a.  Enter  the  graph  on  the  firebrand  initial  height  scale  [from  step  4),  drawing 
a  vertical  line  up  to  the  curve  labeled  with  the  treetop  height  recorded  in  step  5. 

b.  From  the  intersection  point  (interpolated  as  necessary)  draw  a  horizontal 
line  into  the  left-hand  panel,  through  to  the  curve  labeled  with  the  windspeed  de- 
termined in  step  5. 

c.  From  the  intersection  point  (interpolated  as  necessary)  draw  a  vertical 
line  down  to  the  left  horizontal  scale.     There  read  off  the  maximum  spot  fire  dis- 
tance and  record  on  the  worksheet. 

The  dashed  lines  on  the  figure  illustrate  these  steps. 

7.  If  the  terrain  over  which  the  firebrand  would  fly  has  a  substantial  varia- 
tion in  elevation,  the  flat-terrain  spot  distance  can  be  corrected  by  the  use  of  the 
graphs  in  figure  9.     First,  determine  which  description  best  fits  the  location  of 
the  firebrand  source:    midslope  on  the  leeward  side  of  a  ridge  (figure  9A) ,  the 
valley  floor  Cfigure  9B) ,  midslope  on  the  windward  side  of  a  ridge  (figure  9C) ,  or 
ridgetop  (figure  9D) .     Next,  record  the  elevation  difference  from  ridgetop  to  valley 
bottom  and  the  distance  from  ridgeline  to  valley  bottom  as  would  be  shown  on  a  map. 
Divide  the  flat  terrain  spot  distance  from  step  6  by  the  ridge-to-valley  distance 
and  record  the  result.     Use  this  quantity  in  the  appropriate  version  of  figure  9 

to  determine  the  ratio  of  spotting  distance  in  sinusoidal  terrain  to  that  over  flat 
terrain.     Multiply  this  number  by  the  flat-terrain  spotting  distance  from  step  6 
and  record  the  result.     This  result  is  the  maximum  spot  fire  distance  as  would  be 
shown  on  a  map . 


^Albini,  F.  A.  and  R.  G.  Baughman.     Estimating  windspeeds  for  predicting  wildland 
fire  behavior.     USDA  For.  Serv.  Res.  Paper  INT-211,  12  p.     Intermt.  For.  and  Range 
Exp.  Stn.,  Ogden,  Utah. 
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WORKSHEET  FOR  ESTIMATING  MAXIMUM  SPOT  FIRE  DISTANCE 
Describe  tree(s)  assumed  to  torch  out. 

Species   DBH  (cm-in)    Height  (m-ft)     Number  simultaneous  

t     (use>  in  step  2)     f  (use  in  step  4)+  (use  in  step  3,  fig  6)'^ 

See  table  II  for  number  of  figure  to  use  in  step  2. 

The  steady  flame  from  one  burning  tree. 

From  fig  (3A-4A-5A) :  Flame  height  (m-ft) .       From  fig  (3B-4B-5B) :  Duration  


3.     From  fig  6:      Flame  height  multiplier  .     From  fig  6:     Flame  duration  multiplier_ 

Adjusted  values:     steady  flame  height  (m-ft).  Steady  flame  duration_ 


(use  in  step  4)  i  (use  in  step  4,  fig  7) + 

4.  Tree  height  (step  1)   ^  Adjusted  steady  flame  height  (step  3)  =   (use  in  fig  7) 

From  fig  7:     Lofted    firebrand    height/steady    flame    height  . 

Multiply  by  adjusted  steady  flame  height  (step  3).     Result:   (m-ft) 

Add  1/2  tree  height  (recorded  in  step  1).     Result:   (m-ft)  (use  in  step  6,  fig  8) 

5.  Windspeed  adjustments.     Record  for  reference: 

At  elevation'   (m-ft)        at  height  (m-ft)      windspeed  is  (km/h-mi/h) 

Avg.  elevation  (m-ft).     Mean  treetop  height  along  firebrand  path  (m-ft) 

r.- ,  ,        ■     r:-  (use  in  fig  2  and  step  6,  fig  8)  t 

Difference  =   (m-ft)   (use  m  fig  2)     ^  ^  r     >      &  ^ 

From  fig  2:    windspeed  relative  to  value  at  same  height  but  at  avg.  elevation  (X) 

Windspeed  reference  height  t  mean  treetop  height  =   (use  in  fig  1  if  greater  than  1) 

See  instructions  for  step  5  (use  of  fig  1  may  not  be  indicated) 

From  fig  1,  or  by  estimate,  windspeed  at  reference  height  relative  to  treetop  (Y) 

Divide  reference  windspeed  by  (X)  and  result  by  (Y)  .     Result  (km/h-mi/h) 

(use  in  fig  8) + 

6.  From  fig  8CA-B)  :     Ma:<:imum  spot  distance  over  flat  terrain  (km-mi) 

7.  Correct  for  terrain  relief  along  flight  path  if  necessary. 
Firebrand  source  is  located  nearest  to: 

Leeward  mi ds lope   valley  floor   windward  midslope   ridgetop  


Use  fig  no:  9A  9B  9C  9D 

Elevation  difference  from  ridge  to  valley  (m-ft) 

Distance  from  ridge  to  valley  as  would  be  shown  on  map  (km-mi) 

Spot  distance  from  step  6      ridge-to-valley  distance  =   (use  in  fig  9) 

From  fig  9(A-B-C-D):     Spot  distance  in  sinusoidal  terrain/flat  terrain  distance  

Multiply  by  maximum  spot  distance  over  flat  terrain  (step  6). 

Result :  km-mi) 

This  is  the  maximum  spot  distance  (as  shown  on  map)  estimate  '^ 


EXHIBIT  1.     SPECIMEN  WORKSHEET 
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Figure  4a. — Height  of  "steady" 
flame  from  huming  of  one  tree 
crown  in  still  air  for  several 
species  common  to  the  North 
Central /Northeast  U.S.A. 
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NUMBER  OF  I DENTICAL  TREES  BURNING  AT  ONCE 


Figiire  6. — Multi-ipVL cation  factors  for  flame  hei-ght  and  dimensionless  dura:t'lcm 
ix)  acGOynt  for  simultaneous  burning  of  several  identical  trees  that  -produce 
a  single  steady  flame  and  buoyant  plume  structure . 
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Figure  7, — Height  of  firebrands  lofted  by  bianing  trees. 
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2.0  1.5  1.0  0.5  0  200         400  600  800  1,000 

MAXIMUM  SPOT  FIRE  MAP  DISTANCE  (MILES)  INITIAL  FIREBRAND  HEIGHT  ABOVE  GROUND  (FEET) 


Figure  8a. — Nomograph  for  predicting  maximim  spot  fire  distance  over  flat  terrain; 
British  Knits  of  measure. 


MAXIMUM  SPOT  FIRE  MAP  DISTANCE  (KILOMETERS)      INITIAL  FIREBRAND  HEIGHT  ABOVE  GROUND  (METERS) 

Figure  8b. — Nomograph  for  predicting  maximim  spot  fire  distance  over  flat  terrain; 
metric  units  of  measure. 
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Figicpe  9a. — Spotting  distance  rmltiplier 
for  sinusoidal  tevvain  relative  to  flat 
terrain.     Spot  source  is  at  midslope  on 
leeward  side.    Reference  wind  profile 
is  the  same  as  the  midslope  profile. 
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Figure  9b. — Spotting  distance  multiplier 
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terrain.    Spot  source  is  at  valley 
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the  same  as  the  midslope  profile. 


I— 

CO 
Q 


< 


<: 
cm 


< 

o 

CO 


o 

Q_ 
CO 


1.4|- 


1.2 


1.0 


0.8 


0.6 


0.4 


^  0.2 


1  THOUS.  FT 


1200  IVl 


ELEVATION  DIFFERENCE 
RIDGE  TO  VALLEY 


WIND 


A 

SOURCE 


I 


J 


0  0.5  1.0  1.5  2.0 

FLAT  TERRAIN  SPOT  DISTANCE  /RIDGE-VALLEY  DISTANCE 


15 


ELEVATION  DIFFERENCE 
RIDGE  TO  VALLEY 


WIND 


A 

SOURCE 


0 


0.5 


1.0 


1.5 


2.0 


FLAT  TERRAIN  SPOT  DISTANCE  /RIDGE-VALLEY  DISTANCE 


Figure  9c, — Spotting  distance  rmltiplieT 
for  sinusoidal  terrain  relative  to  flat 
terrain.     Spot  source  is  at  midslope  on 
windward  side.    Reference  wind  profile 
is  the  same  as  the  midslope  profile. 
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Figure  9d. — Spotting  distance  multiplier 
for  sinusoidal  terrain  relative  to  flat 
2.0     terrain.     Spot  source  is  at  ridgetop. 
Reference  wind  profile  is  the  same  as 
the  midslope  profile. 
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EXAMPLES 


Two  examples  are  given  to  help  familiarize  the  reader  with  computation  proce- 
dures.    From  the  statements  of  the  problems  one  can  obtain  all  the  information  neces- 
sary to  complete  the  worksheets  as  shown.    The  reader  is  encouraged  to  reconstruct 
the  process  to  see  how  each  number  is  used. 

Example  1,     A  recurrent  fire  control  problem  on  the  Coconino  National  Forest 
near  Flagstaff,  Arizona,  is  spotting  when  thickets  of  small  Ponderosa  pine  torch  out 
under  strong  wind.     Estimate  maximum  spot  fire  distance  if  a  group  of  20  trees,  each 
6  inches  in  diameter,  torches  out  when  the  windspeed  is  20  mi/h,  measured  at  the 
nearest  Fire-Danger  Rating  station  (i.e.,  a  20-foot  tower  anemometer  in  a  clearing). 
The  trees  that  torch  out  are  about  50  feet  tall;  the  terrain  can  be  considered  flat 
with  broken  cover  of  50-foot  trees. 

This  problem  is  shown  in  the  completed  worksheet  of  exhibit  2.     Note  that  in 
step  5  no  windspeed  correction  for  height  above  terrain  was  necessary,  since  the 
20-foot  tower  windspeed  is  close  to  the  "mean  treetop  height"  (25  feet) . 

WORKSHEET  FOR  ESTIMATING  MAXIMUM  SPOT  FIRE  DISTANCE 

1.  Describe  trees (s)  assumed  to  torch  out. 

Species  Pn^e/^mta   PiKe         DBH    (p    (cm  Height  SO    (m-(^t})  Number  simultaneous  ZQ 

t  (use  in  step  2)  t  (use  in  step  4)+  (use  in  step  3,  fig  6)1- 

See  table  II  for  number  of  figure  to  use  in  step  2. 

2.  The  steady  flame  from  one  burning  tree. 

From  fig  ((3a)4A-5A)  :  Flame  height    2>9     (m^^tj)       From  fig  (3b)-4B-5B)  :      Duration  8>2. 

From  fig  6:       Flame  height  multiplier  3.2.^"  ■     From  fig  6:     Flame  duration  multiplier   ,  5"S" 

3.  Adjusted  values:     Steady  flame  height  (m-^^.  Steady  flame  duration   ^>  S" 

(use  in  step  4)   t  (use  in  step  4,  fig  7)+ 

4.  Tree  height  (step  1)       Adjusted  steady  flame  height  (step  3)  =     .5-3   _(use  in  fig  7) 
From  fig  7:     Lofted  firebrand  height/steady  flame  height  6   . 

Multiply  by  adjusted  steady  flame  height  (step  3).     Result:  6  2.0  (m-^^ 

Add  1/2  tree  height  (recorded  in  step  1).    Result:        (^5^      (m-^t^(use  in  step  6,  fig  8) 

5.  Windspeed  adjustments.     Record  for  reference: 

At  elevation  —  (m-ft)        at  height   3.Q  (m-^j)      windspeed  is  vLO  (km/h^mi/^ 

Avg.  elevation  ^  (m-ft).     Mean  treetop  height  along  firebrand  path     2.5"  ('"'^^ 

n- r  r  .c-     ',^     (use  in  fig  2  and  step  6,  fig  8)  t 

Difference  =  —  (m-ft)   (use  m  fig  2)     ^  ^  f     >      e,  j 

From  fig  2:    windspeed  relative  to  value  at  same  height  but  at  avg.  elevation        1.  (X) 

Windspeed  reference  height  t  mean  treetop  height  =  _0^_8__(use  in  fig  1  if  greater  than  1) 

See  instructions  for  step  5  (use  of  fig  1  may  not  be  indicated) 

From  fig  1,  or  by  estimate,  windspeed  at  reference  height  relative  to  treetop  j[  (Y) 

Divide  reference  windspeed  by  (X)  and  result  by  (Y)  .     Result        /tO  (km/h-^i/R\ 

(use  in  fig  8)+   

6.  From  fig  8^A^B) :     Maximum  spot  distance  over  flat  terrain   ,  (km^mj) 

7.  Correct  for  terrain  relief  along  flight  path  if  necessary. 
Firebrand  source  is  located  nearest  to: 

Leeward  midslope   valley  floor   windward  midslope   ridgetop_ 


Use  fig  no:  9A  9B  9C  9D 

Elevation  difference  from  ridge  to  valley  (m-ft) 

Distance  from  ridge  to  valley  as  would  be  shown  on  map  (km-mi) 

Spot  distance  from  step  6  t  ridge-to-valley  distance  =   (use  in  fig  9) 

From  fig  9(A-B-C-D):    Spot  distance  in  sinusoidal  terrain/flat  terrain  distance  

Multiply  by  maximum  spot  distance  over  flat  terrain  (step  6) . 

Result ;  ■  ^  ^  (km 

This  is  the  maximum  spot  distance  (as  shown  on  map)  estimatet 
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EXHIBIT  2.     EXAMPLE  1  COMPLETED 


Example  2.     While  a  wildfire  is  topping  out  on  a  ridge  3,000  feet  above  a  river 
valley,  fire  scouts  note  that  groups  of  5  to  10  large  Douglas-firs  are  torching  out 
sporadically.     These  trees  are  near  the  ridgetop  and  flying  embers  may  be  borne 
across  the  river,  about  one-half  mile  from  the  ridge.    The  trees  are  probably  100 
feet  tall,  and  appear  to  be  in  the  24-30  inch  d.b.h.  range.     The  fire  weather  fore- 
caster estimates  the  windspeed  over  the  ridgetop  to  be  15-20  mi/h  (presumably  at 
treetop  height) .     The  hillside  is  uniformly  forested  with  trees  averaging  75  feet 
in  height.     As  fire  behavior  officer,  you  are  asked  if  you  think  a  crew  or  two  should 
be  dispatched  to  the  far  side  of  the  river  to  extinguish  spot  fires.     Assuming  that 
the  humidity  is  low  so  potential  spot  fires  can  be  easily  ignited,  what  would  your 
answer  be? 

This  example  is  shown  in  exhibit  3.     Note  that  the  "high  end"  numbers  are  used 
in  each  case  to  get  a  maximum  estimate  of  spot  distance.     So  the  trees  are  assumed 
to  be  30  inches  d.b.h.,  and  groups  of  10  are  assumed  to  bum  at  one  time. 

WORKSHEET  FOR  ESTIMATING  MAXIMUM  SPOT  FIRE  DISTANCE 

1.  Describe  treeCs)  assumed  to  torch  out. 

Species  OOOtllaS''  f^fY  DBH  30  [cm^i^    Height  ^QQ  (m-^^   Number  simultaneous  /Q 

»    1-  (use  in  step  2)  +  (use  in  step  4)t  (use  in  step  3,  fig  6)  + 

See  table  II  for  number  of  figure  to  use  in  step  2. 

2.  The  steady  flame  from  one  burning  tree. 

From  fig  ^3^4A-5A)  :  Flame  height   ^7    Cm-^t^       From  fig  ^-A^-S^):      Duration  V.Q 

3.  From  fig  6:      Flame  height  multiplier  3  » S"  •     From  fig  6;     Flame  duration  multiplier  ^03 

Adjusted  values:    steady  flame  height  2.jL^  (J"-^^-  Steady  flame  duration jj  »  S" 

(use  in  step  4)  t  (use  in  step  4,  fig  7)+ 

4.  Tree  height  (step  1)  ^  Adjusted  steady  flame  height  (step  3)  =    « (use  in  fig  7) 
From  fig  7:     Lofted  firebrand  height/steady  flame  height  ^»  ^ 

Multiply  by  adjusted  steady  flame  height  (step  3).     Result:   (m-(f^^ 

Add  1/2  tree  height  (recorded  in  step  1).     Result:  (m-^£)  (use  in  step  6,  fig  8) 

5.  Windspeed  adjustments.     Record  for  reference; 

At  elevation        3000      ('"-^D^       at  height  7S"  ('"-^j^      windspeed  is     XO  (km/h(mi/]i^ 

Avg.  elevation      f  S"0O       (''"-^Zj  ■     Mean  treetop  height  along  firebrand  path      7^"  1^-^^) 

_  (use  in  fig  2  and  step  6,  fig  8)  t 

Difference  =         I SOQ      (m-@)   (use  in  fig  2) 

From  fig  2:    windspeed  relative  to  value  at  same  height  but  at  avg.  elevation       f  *  ^  (X) 
Windspeed  reference  height  t  mean  treetop  height  =     f.Q      (use  in  fig  1  if  greater  than  1) 
See  instructions  for  step  5  (use  of  fig  1  may  not  be  indicated) 

From  fig  1,  or  by  estimate,  windspeed  at  reference  height  relative  to  treetop  j[  (Y) 

Divide  reference  windspeed  by  (X)  and  result  by  (Y)  .     Result       f  ^  (km/h-^^Ii/hS 

(use  in  fig  8)  -t  \._-^ 

6.  From  fig  8^^B)  :    Maximum  spot  distance  over  flat  terrain        O'  (p  ^  (km^i)^ 

7.  Correct  for  terrain  relief  along  flight  path  if  necessary. 
Firebrand  source  is  located  nearest  to: 

Leeward  midslope   valley  floor   windward  midslope   ridgetop 

Use  fig  no:  9A  9B  9C  /"9D 

Elevation  difference  from  ridge  to  valley     3O0O  (m-^^ 
Distance  from  ridge  to  valley  as  would  be  shown  on  map     Q » (km-^T)) 

Spot  distance  from  step  6  i  ridge-to-valley  distance  =       I,  3  (use  in  fig  9) 

From  fig  9(A-B-C^^:    Spot  distance  in  sinusoidal  terrain/flat  terrain  distance  

Multiply  by  maximum  spot  distance  over  flat  terrain  (step  6) . 

0.(c1 


Result:  •  ^  I   (kmfni^ 

This  is  the  maximum  spot  distance  (as  shown  on  map)  estimate  + 

EXHIBIT  3.     EXAMPLE  2  COMPLETED 
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SYNOPSIS  OF  SUBMODEL  CONTENT 


The  individual  submodels  of  the  prediction  scheme  are  described  in  detail  in 
the  series  of  appendixes  following.     This  section  describes  how  the  submodel  results 
are  fitted  together  and  discusses  some  key  assumptions  used  in  developing  the  pre- 
diction scheme, 

The  steady  flame  structure  submodel   (appendix  A)  is  a  crude  phenomenological 
description,  predicting  flame  height  to  be  proportional  to  the  rate  of  fuel  consump- 
tion raised  to  the  power  2/5.     Using  this  model  to  describe  the  flame  from  a  burning 
tree  requires  that  further  assumptions  be  made.     The  assumptions  used  here  are: 

1.     The  fuel  consumption  rate  is  approximated  as  70  percent  of  the  dry  weight 
of  the  crown  foliage  divided  by  a  burning  time  fixed  by  the  surface/volume  ratio  of 
the  foliage.     The  factor  0.7  is  an  arbitrary  estimate  of  the  amount  of  foliage 
burned  during  the  "steady''  flaming  period  when  the  entire  tree  crown  is  involved  in 
fire.     The  burning  time  formula  is  based  on  an  em.pirical  rule  determined  by  Anderson 
C1969).     Crown  foliage  weight  can  be  predicted  for  many  tree  species  from  diameter 
at  breast  height   (Keays  1971;  Brown  and  others  1977). 

2.  Tlie  height  of  tne  uase  of  the  flame  is  not  well  defined.     Tlie  steady  flame 
model  pictures  tne  flow  of  combustible  fuel  gas  through  a  iiorizontal  orifice.  The 
height  of  such  an  equivalent  fuel  supply  plane  is  clearly  between  the  ground  and  the 
top  of  tne  tree.     Tnis  height  represents  the  vertical  displacement  of  the  origin  of 
the  coordinate  system  in  which  the  flame  structure  is  descrioed.     In  this  presenta- 
tion the  origin  is  taken  to  be  at  half  the  tree  height. 

3.  Mien  several  trees  are  consumed  nearly  simultaneously  to  produce  a  single 
large  flame  structure,  the  firebrand- lifting  capability  is  increased.     Tlie  flame 
height  increase  factor  is  taken  to  oe  N^'^^,  implying  merely  a  pooling  of  the  com- 
bustible gases  from  N  burning  trees.     Flames  from  neighboring  fires  will  merge  into 

a  single  structure  when  the  flame  heights  are  significantly  greater  than  the  distance 
separating  them  (Tliomas  and  others  1965).     By  assuming  that  the  flames  always  merge, 
we  clearly  follow  the  intention  of  analyzing  the  "worst  case."      The  merging  of  the 
flames  results  in  a  reduction  of  the  dimens-ionless  flame  duration,  because  the  nor- 
malizing time  is  the  time  required  for  gas  to  pass  through  the  flame  structure.  As 
the  flame  height,  z„,  increases,  this  travel  time  increases  as  z„-^/^.     But  the  ac- 

r  r 

tual  burning  period  stays  constant,  so  the  dimensionless  measure  of  flame  duration 
decreases  as  N"-"-/^. 

The  buoyant  plume  above  the  steady  flame  is  modeled  (appendix  B)  as  a  self- 
similar,  radially  symmetrical,  turbulent  ouoyant  jet.     The  velocity  and  buoyancy 
flux  are  matched  at  the  flame  tip.     Tne  rest  of  the  flow  description  is  obtained 
immediately  by  application  of  well-established  scaling  laws.     So  closely  related 
are  tne  flame  structure  and  plume  structure  models  that  the  flame  height  serves 
to  scale  the  plume  geometry.     The  characteristic  gas  transit  time  also  serves 
naturally  to  nondimensionalize  the  time  history  of  tne  collapse  of  the  plume  when 
the  fire  bums  out.     So  if  the  flame  structure  model  is  revised,  the  plume  flow 
model  will  be  automatically  adjusted  if  the  flame  tip  gas  velocity  scales  as  z 

r 

The  firebrand  burning  rate  model  (appendix  C)  used  in  this  theory  is  quite 
similar  in  form  to  one  inferred  uy  others  (Lee  and  Hellman  1970  ,  Tarifa  and  others 
1965,  1967).     Numerical  comparisons  are  difficult,  however,  and  there  is  a  great 
deal  of  scatter  in  the  supporting  data.     This  model  predicts  that  a  freely  falling, 
burning  cylinder  or  plate  will  fall  at  a  velocity  (relative  to  the  vertical  com- 
ponent of  windspeed)  that  declines  linearly  in  time.     This  prediction,  coupled  with 
the  results  from  a  wind  field  model  (appendix  E)  that  streamlines  parallel  the 
terrain,  leads  to  a  great  simplification  in  determining  the  maximum  possible  spot 
fire  distance. 
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While  firebrand  materials  other  than  cylindrical  stem  sections  are  obvious  can- 
didates for  analysis   (Clements  1977) ,  only  minor  differences  in  spot  distance  would 
result  unless  the  burning  rate  behavior  departs  substantially  from  that  used  in  this 
model.     The  height  to  which  a  particle  will  be  lofted  by  the  flame/plume  flow  (ap- 
pendix D)  depends  upon  the  parameter  group  pD/C^,  where  p  is  the  particle  mean  mass 

density,  D  is  a  linear  dimension  of  the  particle  parallel  to  its  velocity  relative 
to  the  gas  flow,  and  Cj-^  is  the  particle's  drag  coefficient.     The  time  history  of 

this  same  parameter  group  determines  how  far  it  will  be  borne  by  a  given  wind  field. 

Stemwood  segments,  needles,  bark  flakes,  and  seed-cone  scales  will  behave  quite 
similarly  in  this  regard  and  a  cylinder  makes  a  good  generic  representative.  Twigs 
with  foliage  attached,  open  whole  seed  cones,  and  clumps  of  moss  are  examples  of 
potential  firebrands  that  might  behave  differently.     But  even  open  seed  cones  would 
have  substantially  higher  initial  falling  velocities  than  small-diameter  stemwood 
segments  (Clements  1977)  so  would  not  be  carried  nearly  so  high  in  the  flow  and  could 
not  travel  as  far.     Moss  clumps,  on  the  other  hand,  could  be  lofted  high  into  the  air 
because  the  mean  density  of  a  highly  porous  clump  is  very  low.     But  the  "staying 
power"  of  a  burning  moss  clump  would  be  determined  by  the  surface/volume  ratio  of 
the  material  and  it  usually  would  be  consumed  before  returning  to  the  ground. 

A  twig  with  foliage  attached  might  be  able  to  outdistance  a  simple  wood  cylin- 
der, however.     During  the  lofting  phase  of  its  flight  it  could  have  a  low  value  of 
pD/Cp  because  of  a  porous  overall  structure.     But  the  foliage  might  soon  be  burned 
away,  leaving  a  wood  cylinder  with  a  substantially  higher  value  of  pD  than  the  com- 
posite structure  had.     This  would  allow  the  particle  sufficient  lifetime  to  return 
to  the  ground  before  being  consumed  and  it  would  start  falling  from  a  height  that 
it  could  not  achieve  without  the  added  drag  area  afforded  by  the  foliage.     Such  a 
"two  stage"  firebrand  presents  a  difficult  analytical  challenge--one  that  has  not 
been  met  fully  here.     A  cursory  attempt  to  establish  a  range  increase  failed. 

Samples  of  live  twigs  with  needles  attached  were  gathered  and  analyzed  by  the 
author.     Four  species  of  tree  were  represented  in  the  samples:     lodgepole  pine, 
ponderosa  pine,  Douglas-fir,  and  Engelmann  spruce.     The  needle  lengths,  number  of 
needles  per  unit  length  of  twig,  twig  diameters,  needle  weights,  and  twig  weights 
were  determined  for  about  10  twig  segments  for  each  species.     The  weight/drag  ratios 
(or  pD/Cp)  of  the  composite  structures  were  calculated  using  average  values  for  the 

variables  needed.     Using  these  data  in  the  particle  height-vs-time  formulas  (appen- 
dix D)  resulted  in  initial  particle  heights  that  were  too  great  to  allow  the  twig 
to  return  to  the  ground  before  burning  out,  for  all  the  examples  tried.     In  other 
words,  no  spot  distance  advantage  was  found  for  the  composite  firebrands.  This 
does  not  mean  that  all  combinations  of  twig  diameter,  needle  length,  nimber  of 
needles  per  unit  of  twig  length,  etc.,  would  not  go  farther  than  a  stem  segment. 
But,  for  the  species  sampled,  the  necessary  combination  was  not  found.     No  further 
effort  was  spent  on  "two  stage"  firebrand  structures. 

The  structure  of  the  wind  field  that  carries  the  burning  brand  has  a  powerful 
influence  on  the  distance  it  travels.     For  flight  over  forested  flat  terrain  a 
simple  logarithmic  variation  of  velocity  with  height  suffices  for  a  windspeed  model. 
Using  the  profile  shown  in  figure  1,  the  path  of  the  burning  brand  is  readily  com- 
puted.    In  this  computation  the  logarithmic  profile  of  windspeed  with  height  is  ex- 
tended below  treetop  height  to  zero  velocity.     A  maximum  distance  would  be  predicted 
by  assuming  the  windspeed  below  treetop  height  to  be  constant.     At  the  other  extreme, 
one  might  consider  the  crown  layer  to  be  impenetrable  by  a  firebrand  and  set  the 
velocity  to  zero.     Neither  extreme  is  realistic,  but  the  extended  logarithmic  pro- 
file interpolates  between  them,  so  is  used  for  convenience. 
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For  wind  over  uneven  terrain--e . g . ,  blowing  across  a  ridge-valley  system--it  is 
necessary  to  model  the  variation  of  windspeed  with  horizontal  distance  (fig.   2)  and 
to  include  the  vertical  wind  component  as  well.     The  crude  model  developed  in  appen- 
dix E  includes  these  variations,  but  omits  some  well  established  features  of  such  a 
wind  field: 

1.  The  variation  of  windspeed  and  direction  with  height  commonly  called  the 
Ekman  spiral   (see,  e.g..  Brown  1974}   is  not  included  in  the  model.     An  isobar  at 
constant  height  is  used  in  the  model  as  a  "lid"  on  the  flow  beneath  it,  as  done  by 
Fosberg  and  others   (1976)   and  friction  is  ignored  down  to  a  thin  layer  close  to  the 
surface.     In  this  surface  friction  layer  the  logarithmic  profile  used  in  the  flat 
terrain  model  is  employed,  with  a  characteristic  horizontal  velocity  that  depends 
on  local  terrain  height. 

2.  The  "rotor"  or  recirculation  flow  often  observed  on  the  lee  side  of  a 
ridge  (Schroeder  and  Buck  1970}  is  not  included  in  the  wind  model.  There  exist 
mathematical  models  of  windflow  over  terrain  features  that  include  "rotor"  and 
"lee  wave"  effects  (Alaka  1960)  but  they  are  complex  and  demand  more  data  than 
just  the  terrain  profile.  Trial  use  of  the  simplified  wind  model  presented  here 
will  reveal  whether  or  not  the  simplification  has  degraded  its  accuracy  beyond 
utility  in  rough  terrain. 

3.  The  friction  length  scale  of  the  flow  in  the  friction  layer  is  treated  as 
a  constant  in  this  model.     Its  value  is  set  to  0.1313  times  the  average  treetop 
height,  based  upon  measured  data  in  relatively  flat  terrain  with  a  uniform  forest 
cover  (Tanner  and  Pelton  1960;  Marunich  1975).     This  may  prove  to  be  a  weak  point 
in  the  velocity  profile  model  because  uniform  forest  cover  over  nonuniform  terrain 
with  significant  elevation  change  is  not  the  usual  case.     The  user  is  asked  to 
select  an  "average"  canopy  top  height  when  assessing  maximum  spotting  distance, 
but  no  guidance  is  given  for  determining  the  appropriate  average  when  forest  cover 
and  open  areas  must  both  be  considered. 

Determination  of  the  trajectory  of  the  firebrand  particle  carried  by  the  model 
flow  fields  is  straightforward.     Although  some  approximations  are  made  in  the  ana- 
lyses, each  can  be  justified  quantitatively.     It  is  felt  that  the  approximations 
and  model  uncertainties  in  the  burning  rate  and  wind  field  descriptions  dominate 
errors  introduced  in  calculating  the  flight  of  the  particle.     The  lofting  phase  of 
the  trajectory  (appendix  E)   is  treated  less  precisely  than  the  transport  phase 
(appendix  F)  but  in  neither  case  should  the  integration  of  the  equations  of  motion 
introduce  error  as  large  as  the  uncertainties  in  the  models  on  which  they  are  based. 
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APPENDIX  A:  FLAME  STRUCTURE  MODEL 


In  this  appendix  a  simple  phenomenological  model  is  presented  for  the  structure 
of  a  large  flame  burning  in  still  air.     Approximate  expressions  are  developed  for 
tlie  mean  length  Q^eight)  of  the  flame,  the  velocity  of  the  gas  flow  within  the 
flame,  and  the  dynamic  pressure  profile.     To  derive  these  expressions  a  model  of 
such  a  flame  is  postulated  and  analyzed.     Although  the  model  is  obviously  over- 
simplified and  incorrect  in  detail,  the  gross  relationships  sought  do  not  justify 
a  more  elaborate  model .     And  the  transparency  of  the  model  allows  one  to  see  the 
influence  of  various  parameters   (and  assumptions)  on  the  quantities  sought. 

I.     Conceptual  Model 

Picture  an  erect  right  circular  cone  as  a  boundary  in  space.     This  cone  marks 
the  region  within  which  exists  gaseous  fuel.*^     The  fuel  flows  into  the  region 
axially  across  the  circular  area  at  the  base  and  is  consumed  at  the  surface  of  the 
cone  as  it  makes  contact  with  the  air.     The  tip  of  the  cone  is  the  tip  of  the  flame. 
The  cone  of  fuel  gas  is  surrounded  by  a  cylindrical  sheath  of  varying  radius  within 
which  flow    tlie  products  of  combustion  plus  entrained  ambient  air.     In  this  annular 
sheath  the  gas  flows  upward,  as  does  the  fuel  gas.     Outside  this  sheath,  ambient 
air  flows  radially  inward  toward  the  vertical  axis  of  the  flame.     The  following 
symbols  will  be  used  to  describe  the  flows  and  boundaries  discussed; 

Symbo Is 

p  =  gas  density  (mass/unit  volume} 
r  =  radial  distance  from  central  axis 
z  =  vertical  distance  from  base  of  flame 
u  =  inward  radial  velocity 
w  =  upward  axial  velocity 

^  =  mean  mass  ratio  of  air  to  fuel  producing  flame-product/air  mix 

n  =  entrainment  velocity  ratio  Cu/w} ,  a  constant  with  height  by  assumption 

g  =  acceleration  of  gravity 

SubsGV^-pts 

f    refers  to  fuel  gas   Ccondition  or  boundary  location} 
a    refers  to  ambient  air  Ccondition  or  boundary  location} 
p    refers  to  products  of  reaction  (plus  admixed  air} 
o    refers  to  conditions  at  the  base  of  the  flame 
F    refers  to  conditions  at  the  tip  of  the  flame 


The  gaseous  fuel  consists  of  the  products  of  volatilization  and  pyrolysis  of 
the  solid  fuel  being  burned,  plus  some  combustion  products  and  some  poorly  mixed 
air.     In  this  simplified  model  the  mixture  is  idealized  as  pure  gaseous  fuel. 
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II.    Ma thema tioal  Mode  I 


Using  these  symbols  we  can  write  the  conservation  equations  of  the  gross  flame 
structure.  The  entrainment  assumption  CTliomas  1963;  Tnomas  and  others  1965;  Turner 
1973)  is  made  here  in  the  form: 

u  =  nw  (Al) 

and  w  is  assumed  to  be  a  function  of  z  only.     The  conservation  of  mass  and  momentum 
are  readily  expressed  in  terms  of  an  intermediate  variable,  m,  representing  the 
mass  in  a  unit  iieight  of  the  entire  flame  structure,  divided  by  tt  ; 

m  =  Pp(r^^  -  r/)  +  p^r/,  (A2) 

The  conservation  of  mass  relationship  is  then: 

^(mw)  =  2p  r  u.  (A3) 
dz  a  a  ^  ■' 

Tlie  air/fuel  ratio  involved  in  forming  the  "product"  gas  (not  chemical  stoichiomet ry) 
is  specified  oy  the  parameter  (j).     Holding  ^  constant  leads  to: 

2  2 

mw  +  (f)p_^wr_p    =  CI  +  <j5)Cp£wr_p  )^.  (A4) 

Tlie  conservation  of  momentum,  including  the  effects  of  buoyancy  and  "drag"  due  to 
air  entrainment  is  expressed  by: 

mw^  =  gCr^  p^  -  m)  -  2p^r^uw.  (A5) 

These  equations  are  essentially  the  same  as  those  used  by  others  (Steward  1964; 
Morton  1965;  Nielsen  and  Tao  1965)  in  models  of  flame  size. 

Using  equations  (A2)  and  (A4)  one  can  derive  the  following  useful  intermediate 
relationship : 

2  _  r^f /1+^x      llmw  ^^r  2. 


wr      =  <^ — ( — — ^)  -  -p-  -   ( — — ^)(  IJCwr-  )   .  (A6) 

^  p       ^  J   f  ^  ^p 

Here  it  has  also  been  assumed  that  the  inner  core  of  fuel  is  at  constant  density: 

Pf  =  (pPo-  (A7) 


Using  equation  (A6)  in  equation  (A5)  multiplied  through  by  w  gives: 

fVPpW    I     ^1    V    (j)  Ap^ 

In  a  shorthand  form  this  can  be  written  as 


2  dw 

mw    -J—  =  P  &i  \ 
dz  a^lp^Vp 


fmw/^f /I  +  l\      /I  +  (l)x/^f      ,\  ,      2.\  ^  2  , 

fcl~(-T^)  -        -  (-T-^)(-  -  ^)  ^^^f  \]  -  -  2p^r^uw.  (A8) 


d  2 

w— (mw  )  =  g(amw  -  3(mw)  )  (A9) 
dz  o 

where  a  and  3  are  constants  given  by: 
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1      '  p       1  p 

a  =  {^^)—  -  -  —  -  1  (AlO) 
9       -p       —  f 

3  =  (^-^)^(—  -  !)■  (All) 

Note  tnat  tne  mass  conservation  and  entrainment  equations  combined  can  be  written  as 

d  _  ^  _  ,..l/2,,.„  2,1/2 

4/2 


,  -fmw)  =  2r,p  r  w  =  2^,c  w    "  (\<t 
dz^^  aa  a  a 


=  2riP  w      J  mw  -  — [mwj 

a        ]      c  p 


a 


r  1^/2 

=  2Ti<P^w(Ca  +  l}mw  -  3(1150^)  f  CA12) 


0 


or 

1/2  d 


m 


^(mO  =  2JCm^0c^((a  +  l}m;s'  -  SCmvOj  |      •  (A13) 


Multiplying  through  oy  [m.w)  makes  the  left  hand  side  a  perfect  differential: 

d         2  ^ 

^((mw  )  )   =  2gmK^^m;N  -  S[m^s•)^).  (A14] 

Equations   (A13}  and  (A14)  are  no\s"  furt'her  simplified  by  noting  that,  over  most 
of  tne  length  of  tne  flame, 

amvs'  >>  s^mu^Qj  CA15) 

an  inequality  tnat  arises  from  tiie  large  ratio  of  air  to  fuel  noted  for  natural 
fires  (Thomas  1963)  .  By  eliminating  the  second  term  from  the  righthand  side  of 
CA13)  and  CA14}  one  arrives  at^ 

m^''^  -^Ci^O  =  2n((a  +  13c^)^/^mw  (A16) 
2 

^(Cmw^)    )  =  2agCmw3^  (A173 

Substituting  the  expressions: 

P  =  miv;        Q  =  m\<^  (AIS) 
into  these  equations  gives: 


^Actually,  it  is  not  necessary  to  make  this  approximation  to  proceed  with  the 
analysis,  but  tne  difference  resulting  is  minor  and  therefore  not  justifiable  in 
view  of  tne  level  of  accuracy  of  this  model. 
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(P/Ql/2^d^^  2n((a  +  1)P^)^/^P  (A19) 
2agP^  (A20) 


dz 

dQ^  _  o.„d2 
dz 


or 


^  =  2n  ((a  +  l)p^j      Q  (A21) 
^  =  agP^/Q.  (A22) 
Dividing  one  equation  by  the  other  gives: 


or 

((a  +  l)p^) 

I.  n((c  +  1)P^)  ' 

III.    Flame  Lengtn 


This  last  expression  can  be  used  directly  in  equation  CA21)  to  provide  a 
formula  for  flame  length.     Note  that 


^  =  2n  (Ca  .  Dp  r    ^    5   og  ^  p3  j 


dz 


1/5 

(A26) 


a> 

Evaluating  tiie  constant  S  at  the  base  of  the  flame  gives 


o     i)/2      -  Cmw) 
c        /-     2,  oag  o  rA0  7^ 

so  that   CA26)  can  be  written  in  terms  of  P  normalized  by  its  value  at  the  flame 

base,  P  .  in  the  form 
o 


d(P/P  J  2n((a  +  l)p  ) 


1/2 


 :  7175  =  — — uY^—  ^^^^ 

m         /  ^     j>        \  1  m 


((a  -  Dp  ^<  ' 
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Now,  since 

F 

tile  flame  length  formula  is  : 

1/2  /-l+c})  1/2 


2n  -^-^^^  )  =    /     {l  ^  12-  ^1/2^^^    -  ^^2°^ 

\         0/  1^  n((a+l)p^  w 


The  second  term  in  the  radical  in  the     integrand    is  much  larger  than  unity  over 
most  of  tne  length  of  the  flame  so  that,  to  a  level  of  approximation  consistent  with 
this  analysis  we  have: 

r     .  1     1/2  /  \  1/2     2  r  l  +  < 


o  agm 


/     r         T.\l/2  2 
12  MPa'^"  "  %   \     5,,,  ,,2/5 


5  1/2  j  .2 

agm 
^  o 


((1  +  ^r""  -  i).  (A31) 


Tliis  result  is  very  similar  to  that  outained  by  Steward  (1964)  and  exhibits  the 
power  law  found  by  other  investigators  (Vvohl  and  others  1949;  Putnam  and  Speich 
1963;  Putnam  1965} . 

Thomas   (1963)  presents  experimental  data  in  dimensionless  form,  with  which  this 
formula  can  be  compared.     In  the  notation  used  here,  his  dimensionless  form  can  be 
written  as 

where  K  and  N  are  empirical  constants.     Noting  that 

%  =  ^Pf^f^^o  =  ^Pf^a^^o^  ^^^^^ 
equation  (A31)  can  be  cast  in  the  form  of  (A32) ,  with  the  result: 

N  =  I  (A34) 

^  4ri(a  +  \y'^ 

The  numerical  value  of  K  for  natural  fires  is  to  be  estimated  by  using  appropriate 
values  for  ^,  n ,  and  the  density  ratios  appearing  in  the  parameter  a.  Typical 
values  might  be  (Thomas  1963) 

P  /p    =  8/3;    p  /p^  =  5/3;  (})  =  30 ;    n  =  0.075.  (A36) 
a    p  at 
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Using  these  values  we  find 

K  =  12.4  (A37) 

which  value  is  aoout  half  of  that  reported  by  Thomas  (1963)  for  the  experiments  of 
Putnam  and  Speich  (1963)  using  city  gas,  adjusted  for  the  ratio  of  heat  content  of 
wood  volatiles  to  that  of  methane.     The  exponent  (N)  of  2/5  is  the  dimensional ly 
correct  value  and  is  a  good  correlation  value  for  the  data  of  Putnam  and  Speich, 
although  Thomas  reports  better  agreement  with  a  variety  of  wood  criu  data  using  an 
exponent  of  0.61.     Note  that  if  one  replaces  the  lower  limit  of  the  integral  in 
equation  (A31}  witii  zero-- amounting  to  complete  neglect  of  the  fuel  gas  flow  at  the 
base  of  tiie  flame  with  respect  to  the   flow  of  product  gas  at  the  tip  --the 
equation  for  tne  flame  length  becomes 

using  the  constants  given  in  CA36)  above.     This  expression  gives  a  better  fit  than 
K  =  12.4  to  the  experimental  data  given  in  Tliomas '  survey,  and  will  be  used  in  this 
paper. 


IV.     Gas  Velocity  and  Dynamlo  Pressure 

The  mean  gas  velocity  within  the  flame  can  be  expressed  in  terms  of  the  distance 
from  the  flame  base,  using  the  expressions  developed  above.     From  (A25)  and  iA27) 
we  can  write  C^/w^)  in  terms  of  iP/P^) ,  since  Q  is  wP : 

3 

f      o  5/2  (mw)  /       3       \ -^2/5 

w/w^  =  Q/w^P  =     (mw  )^      .  ^  f  _  (_)  _  1        /  (w^P)  .  (A39) 

L  (p^(a  +  1))  '     \    o  /J  / 

Again  neglecting  tne  smaller  terms  in  the  radical  this  can  be  approximated  by: 

^-^4^-,  °—TTn\       ^       •  CA40) 


.12  n    ((„  ,  13     )  1/2 


o 


The  dependence  of  (P/P  )  on  z  is  readily  obtained  from  (A31)  by  replacing  z    by  z 

O  r 

1/2 


and  (1+qj)  by  (P/P^)  •     So  doing  leads  to  the  form 


This  expression  simplifies  also,  since  we  are  concerned  only  with  the  upper 
part  of  the  flame,  where  the  second  term  in  the  radical  on  the  righthand  side  is 
much  greater  than  unity.     Neglecting  the  one  leads  to  considerable  cancellation 
of  factors: 

w  =  (agz/3)^^^.  (A42) 
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The  value  of  a,  from  (AlO) ,  is  very  nearly 

a  =  p^/pp  -  1  (A433 
since  (j)  is  about  30.     This  leads  to  the  form: 

w  =  0.408(2g(p^/Pp  -  l)z)^/^.  (A44) 

The  only  difference  between  this  expression  and  the  empirical  equation  given 
in  Thomas  (1963)  ,  based  on  experimental  work  by  Rasbash  and  others   (1956)  with 
liquid  fuel  fires,  is  the  numerical  factor.     Thomas  reports  a  coefficient  of  0.36, 
or  12  percent  less  than  the  value  here.     The  difference  is  no  doubt  due  to  the  gross 
nature  of  this  model. 

The  mean  density  within  the  flame  structure  is  needed  to  describe  the  dynamic 
pressure  profile.  Average  density  (p)  is  readily  derivable  from  the  variables  al- 
ready used;  since  by  definition  of  average 

-  2  2  2 

p  =  m/r    =  raw/ (wr  )  =  P/ (wr  ).  (A45) 

d.  H.  3. 

From  (A2) ,   (A6) ,   (A7) ,  and  (AlO)  we  have 

wr^  =  ±±±  P  -  3  ^  (A46) 
a         p  p  ^  ^ 

a  a 

so  that 

p/p^  =  P((a  +  1)P  -  aP^)'"^  =    ((a  +  1)   -  3Pq/p)'^  (A47) 

This  expression  simplifies  by  referring  the  mean  density  to  the  value  at  the 
tip  of  tne  flame  and  using  the  scaling  law  for  the  parameter  P  referred  to  its 
value  at  the  flame  tip: 

P/Pp  =   (z/Zp)^/^      P  =   (1  +  (^)P^(z/Zp)^/^  (A48) 

Using  this  form  and  noting  that  p  at  the  flame  tip  is  Pp,  equation  (A47)  can  be 
approximated  for  the  upper  part  of  the  flame  by: 

p/Pp  =  (l  +   (1  -  Pp/p^)(l  -  (Zp/z)^/2)/(}))"^  (A49) 

From  this  form  it  is  apparent  that,  over  most  of  the  upper  part  of  the  flame,  the 
density  can  be  considered  constant  and  equal  to  the  flame  tip  value. 

The  dynamic  pressure,  referred  to  its  value  at  the  tip  of  the  flame  is  quickly 
de  ri  ve  d  n  ow  : 

5 /2 

-V/   =  ^^^K:^)  -))■ 

r     r  r 
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V.     Formulae  Used  in  Firehrand  Lofting  Analysis 


Tiie  formulae  derived  above  can  be  cast  in  practical  form  for  use  in  the  firebrand 
lofting  computations.     Equation  (A38)  gives  the  flame  height,  but  in  terms  difficult 
to  apply.     It  can  be  simplified  if  one  observes  that  the  rate  of  fuel  consumption 
(mass  per  unit  time,  M)  can  be  written  as: 

M  =  up^(r^w)   .  (A51) 
f^  a  o 

Inserting  this  into  (A38)  and  canceling  out  (r  )  gives 

3.  O 

1/2  •  2/5 

Assuming  approximately  standard  atmospheric  conditions 

p     =0.075  Ib/ft^  =1.2  kg/m^  CA53) 


a 


32  ft/s^  =9.8  m/s^  (A54) 


so 


Zp  =  25 (M)^/^  (A55) 

where  the  flame  height  is  in  feet  and  the  fuel  consumption  rate  is  in  pounds  mass 
per  second,  or 

Zp  =  10. 4  (M)^'^^  (A56) 

if  the  units  are  meters  and  kilograms  per  second.     Figure  A-1  compares  predicted 
flame  heights,  using  this  equation,  to  values  observed  in  the  burning  of  wood  cribs. 
The  data  were  transcribed  graphically  from  figures  presented  by  Thomas  (1963).  The 
agreement  between  the  predicted  heights  and  observed  heights  is  quite  satisfactory, 
but  application  of  tne  equation  to  burning  trees  involves  extrapolation  beyond  the 
range  shown  in  figure  A-1. 

For  velocity  all  one  needs  to  do  is  express  the  value  at  height  z  as  a  function 
of  the  flame  tip  value 

w/Wp  =  (z/Zp)^/^  (A57) 

and  calculate  w    from  (A44} .     Here  we  assume  that  the  density  ratio  p  /Pp  is  8/3, 
as  before,^  to  obtain 


°In  other  words,  the  flame  edge  is  defined  as  the  500  C  outline,  about  the 
limit  of  carbon  particle  visiuility,  and  the  product  gas  molecular  weight  is  essen- 
tially tnat  of  air. 


34 


4.2z 


2.  3z 


1/2 
F 

1/2 


£t/s  ,     z „  in  ft 
r 


m/s ,      z    in  m. 

r 


(A5  8) 


The  peak  dynamic  pressure  in  the  flame  occurs  at  the  flame  tip  in  this  model. 


and  the  value  would  be  given  by  q  ,  where 

r 

^F  =    I'A  ~-  ^^a^^F^^  -  Pp/^a^ 


(A59) 


or 


^F  = 


0.00  78  z„      Ibf/ft   ,      z^  in  ft 

r  r 


1.23  z, 


N/m^ , 


(A60} 


z„  m  m. 

r 


For  all  practical  purposes  the  dynamic  pressure  can  be  considered  to  increase 
linearly  with  height. 


Figure  A-1. — Comparison  of 
data  on  flame  heights 
(Laboratory  Wood  Crib 
Burns  J  Thomas  1963) 
with  predictions  from 
present  theory. 
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APPENDIX  B:  THE  BUOYANT  PLUME  ABOVE  A 

STEADY  FLAME 


The  vertical  transport  of  a  firebrand  into  the  prevailing  windstream  may  in- 
clude a  period  of  flight  in  the  buoyant  plume  flow  above  the  flame  from  which  it 
came.  Tliis  "lofting"  phase  can  be  approximated  by  an  upper-limit  calculation,  as 
discussed  in  tiie  text.  To  perform  the  trajectory  calculations  it  is  necessary  to 
describe  the  fluid  flow  field  in  the  buoyant  plume.  From  the  tip  of  the  (transient) 
flame  upward  we  use  iiere  the  scaling  laws  for  a  fully-developed  turbulent  buoyant 
plume  in  still  air. 

The  decay  of  velocity  in  a  buoyant  plume  occurs  through  accretion  of  ambient 
air  into  the  flow.     In  still  air  this  process  is  usually  represented  by  an  "en- 
trainment  parameter"  as  used  in  appendix  A,  giving  the  ratio  of  inflow  velocity 
to  vertical  velocity.     Tlie  entrainraent  parameter  is  of  the  order  of  0.1,  so  use 
of  a  "still  air"  model  clearly  gives  an  upper  limit  to  the  firebrand-lofting  capa- 
bility of  the  plume.     A  plume  in  a  shearing  wind  environment  would  be  eroded  more 
rapidly  witii  height. 

Tlie  fact  that  the  plume  is  treated  as  fully  developed  (thsX  is,  steady  rather 
than  ciianging  witii  time)  is  probably  not  as  conservative  as  it  may  seem.     This  is 
so  because  prior  to  the  "torching"  of  the  tree  or  trees  that  produce  the  plume,  a 
substantial  fire  must  iiave  burned  in  the  same  place  for  a  considerable  period  of 
time.     This  fire  will  have  established  a  convection  plume  that  is  then  merely  am- 
plified by  the  heat  pulse  from  the  burning  tree  crownCs). 

Well-accepted  analysis  of  a  free  buoyant  plume  CTumer  1973)  provides  the 
scaling  laws  needed.     The  simplest  applicable  model  is  one  that  describes  the 
self-similar  flow  field  for  constant  buoyant  flux.  Let 

b  =  mean  radius  of  plume 
G  =  buoyancy  =  gQ  -  P/P  ) 

Si 

w  =  mean  vertical  velocity 

-  2 

F  =  buoyant  flux  =  wb  G  =  constant 


wnere 


g 
P 


acceleration  of  gravity 
mean  density  in  the  plume 
ambient  air  density. 


The  scaling  laws  of  interest  are   (Turner  19  73) 


b 


6 


(Bl) 


w 


1/3  -1/3 


(B2) 


G 


) 


-1/3  -b/3 


(B3) 


where  z  is  height  and  n  is  the  entrainment  parameter. 
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Forming  ratios  of  these  variables  with  their  values  at  the  tip  of  the  flame 

(where  z  =  z  )  we  have: 
r 

w/Wp  =  (Zp/z)^^^  (B4) 

G/Gp  =   (1  -  p/p^)/(l  -  Pp/p^)  =  Czp/z)^/^  (B5) 

The  mass  density  at  the  tip  of  the  flame  is  Pp,  as  in  appendix  A.  The  density  pro- 
file is  given  by 

p/p^  =  1  -   (1  -  Pp/p^)CG/Gp)  =  1  -  (1  -  Pp/p^)CZp/z)^/^  (B6) 
whicli  leads  directly  to  the  dynamic  pressure  decay  law: 

q/qp  =    (p/Pp)(w/Wp)^  (B7) 

q/qp  =  (l  -  (1  -  Pp/p^)(Zp/z)^/^)(p^/Pp)(Zp/Z)2/^  (B8) 
Using  tlie  nominal  density  ratio 

Pp/P^  =  3/8  (B9) 
that  was  used  in  appendix  A,  equation  (B8)  becomes 

q/qp  =   C8/3)(l  -   C5/83(Zp/z)^/^)CZp/z)2/^  (BIO) 


Equation  (BIO)  exhibits  a  maximum,  due  to  the  fact  that  the  density  increase 
rate  exceeds  the  velocity  decrease  rate  in  the  lower  part  of  the  plume.     This  fea- 
ture assures  that  a  firebrand  particle  lifted  within  the  flame  should  be  projected 
above  the  height  of  maximum  dynamic  pressure,  v\fhich  is  at  about  1.6  times  flame 
height.     The  region  oetween  the  flame  tip  and  somewhere  near  this  height  can  be 
thought  of  as  a  transition  zone  between  the  flame  and  the  dissipating  plume  above 
it. 


In  order  to  express  analytically  the  flight  of  a  particle  borne  in  the  flame/ 
plume  flow  structure  it  is  necessary  to  approximate  the  power- law  relationships 
given  above  by  simpler  forms.     The  "transition  zone"  concept  is  valuable  in  for- 
mulating the  needed  approximations.     As  noted  in  appendix  A,  the  mean  fluid  density 
is  nearly  constant  in  the  upper  half  of  the  flame,  while  the  velocity  increases  as 
the  square  root  of  height.     In  the  transition  zone  the  density  increases  substan- 
tially while  the  velocity  changes  very  little.     Above  the  transition  zone,  the 
density  again  rises  slowly  while  the  velocity  gradually  decays.     Simplified  approx- 
imations,    based  on  fitting  the  curve  forms  given  in  this  and  the  previous  appendix, 
are  listed  below : 
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p/p 


(5/8)  (z/Zp3  ,     Zp  <  z  <  1.4  Zp 


(Bll) 


w/w  = 


0.525       ,     z  >  1 .4  z, 


(z/Zp) 


1/2 


.    z  <  z, 


,     Zp  <  z  <  1 . 4  Zp 


(B123 


C5. 963/ (4. 565  +  z/z^))^^^,     z  >1.4  z„, 

r  r 


These  simpler  forms  are  used  in  the  analysis  of  trajectories  of  potential 
firebrands.     The  density  is  underestimated  by  the  approximate  form  and  the  velocity 
is  overestimated,  but  the  dynamic  pressure  is  accurately  described  (fig.  B-1). 
Overest imation  of  the  velocity  is  consistent  with  the  intent  here  to  estimate  an 
upper  bound  for  spotting  distance;  the  misestimation  takes  the  form  of  a  shift  up- 
ward (of  about  2  flame  heights)  of  the  velocity-vs -height  curve. 


1.4  I- 


< 
> 


LU 
> 


LU 

ZD 
CO 

LU 

Ql. 
O 


POINTS  CALCULATED  FROM  EQ.  B  8 
I  (BUOYANT  PLUME  FORMULATION) 


2  3  4  5 

HEIGHT,  MEASURED  IN  FLAME  HEIGHTS 


Figvre  B-1. — Dynamic  pressure  profile  in  and  above  steady  flame ^  showing  approximation 
used  in  calculatinc  firebrand,  liftino  ahilitu. 
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APPENDIX  C:  FIREBRAND  BURNING  RATE  MODEL 


A  firebrand  burning  in  flight  changes  aerodynamic  properties  in  at  least  two 
ways  simultaneously.     It  loses  mass  and  it  loses  volume.     Woody  material  is  pyro- 
lyzed  by  heat  supplied  by  combustion  on  the  outer  surface  of  tlie  piece,  and  is  even- 
tually reduced  totally  to  char.     The  density  of  char  is  less  than  that  of  sound  or 
partially  decayed  wood,  so  a  firebrand  also  loses  density.     By  losing  both  density 
and  volume,  the  firebrand  reduces  its  terminal  fall  velocity,  since  that  velocity 
is  proportional  to  the  square  root  of  the  product  of  density  aiid  linear  dimension 
in  the  direction  of  relative  wind  (Tarifa  and  others  1965} . 

We  are  dealing  with  relatively  small,  light  objects  as  firebrands  here  so  ter- 
minal fall  velocity  is  achieved  in  a  short  distance.     The  subsequent  flight  of  the 
object  is  dictated  by  the  changes  in  terminal  velocity  over  time  (or  distance).  Thus 
the  simplest  model  one  might  construct  to  represent  the  burning  of  firebrands  in 
flight  would  be  one  that  allowed  the  prediction  of  terminal  fall  velocity  directly. 
This  fact  was  recognized  by  Tarifa  and  others  (1967).     The  final  product  of  that 
investigation  was  a  correlation  of  terminal  velocity  with  a  complicated  variable 
that  is  linear  in  time.     The  correlation  curve  is  roughly  bell-shaped. 

The  data  reported  in  these  earlier  investigations  dealt  with  various  object 
shapes  (plates,  cylinders,  and  spheres)  of  differing  sizes  and  made  from  different 
woods.     These  objects  do  not  necessarily  typify  naturally  occurring  potential  fire- 
brands, but  the  fact  that  they  all  yielded  a  similar  behavior  pattern  encourages 
the  use  of  a  "universal"  burning  rate  law. 

Data  taken  at  the  Northern  Forest  Fire  Laboratory'^  on  limbwood  sections  with 
bark  cover  were  used  to  develop  the  model  used  in  this  analysis.    The  burning  rate 
law  is  postulated  from  dimensional  considerations  and  can  be  derived  as  follows. 
Let 

M  =  mass  of  burning  particle 
=  mean  density  of  particle 

D  =  dimension  of  particle  in  direction  of  relative  wind 

(=  thickness  of  plate,  diameter  of  cylinder  or  sphere) 
A  =  cross-sectional  area  presented  to  wind  by  particle 
U  =  velocity  of  wind  over  particle 
p    =  density  of  air. 


Data  reported  by  Muraszew  and  others  (1975) . 
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Tlie  rate  of  mass  loss  is  supposed  to  be  proportional  to  the  rate  of  supply  of 
air  to  tue  surface  (Lee  aiid  hellmaii  1970), 


^  =  -K^p  AU  (CI) 
dt  la 

where        is  a  dimension less  constant.  Since 

M  =  Cp^AD  (C2) 

wnere  C  is  a  geometrical  constant  fixed  by  tne  shape  of  the  particle,  equation  (CI) 
can  oe  written  as 

Tiie  data  on  limuwood  sections   (taule  C-1)  showed  tne  fractional  decrease  in  A 
usually  to  oe  mucii  less  than  tne  fractional  decrease  in  M  over  most  of  the  life- 
time of  ouming  urands,^  so  to  a  fair  approximation  one  might  neglect  the  second 
term  on  the  left-hand  side  of  (C3)  and  obtain 

l-tp^D)  ^        P,U  =  -KP^U.  (C4) 


For  constant  windspeed,  such  as  in  tiie  wind  tunnel  tests,  equation  (C4)  can 
oe  integrated  to  give: 

CP3D)/(P3D)^  =  1  -  Kp^Ut/CP3D)^  (C5) 

wnere  tne  subscript  o  indicates  the  initial  value  of  the  parameter.  This  is  tiie 
form  that  was  used  to  estimate  a  value  for  the  constant  K.  A  constrained  linear 
regression  of  the  form 

y  =  Dx  (C6) 

witu 

y  =  1  -   (P3D)/(p^D)^  (C7) 

X  =  p^Ut/(p^D)^  (C8) 

was  performed  on  tne  data  from  33  tests  on  ovendried  limbwood  sections.  The  variable 
y  was  established  oy  using 

y  =  1  -    ((p^M)/(p^M)J^/2  (C9) 

since  botn  mass  and  density  were  determined  after  fixed  periods  of  burning,  but 
diameter  was  not  directly  measured  except  for  a  limited  series  of  samples. 


8p 

or  cylinders  in  cross  flow,  A  is  tue  product  of  length  and  diameter  (D) . 
While  D  changed,  the  lengtn  was  essentially  constant  uiiless  the  piece  fractured. 
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Table  C-1 . --Firebrand  burning  rate  data.     Data  are  for  ovendry 
limbwood  sanpteSj  each  initially  5  inches  long 
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pD/(pD)     =  density  ratio  times  diameter  ratio. 
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Substantial  scatter  was  evident  in  the  data,  and  two  sources  of  such  scatter 
can  be  identified.     In  some  trials  not  all  of  the  sample  burned.     Wien  only  part 
of  the  lengtii  of  tne  section  burned,  the  final  density  and  mass  would  be  much 
greater  than  in  those  cases  where  the  whole  section  burned.     Also,  in  some  trials 
(especially  for  longer  oum  times)  the  sample  would  bum  through  and  pieces  would 
be  blown  away.     In  tnese  cases  the  final  mass   (but  not  final  density)  would  oe  sig- 
nificantly lower  tiian  in  cases  that  had  wliole  samples  at  the  end  of  the  run.  Hie 
records  showing  wnich  tests  suffered  which  defect  (if  any)  were  not  available^  at 
the  time  of  tnis  analysis,  so  all  data  were  used,  with  the  foreknowledge  that  the 
results  would  include  significant  variation.     Figure  C-1  illustrates  the  scatter 
mentioned. 


Altnough  not  enough  data  are  available  yet  to  quantify  the  process,  the  burn- 
ing of  wood  pieces  under  forced  convection  appears  to  proceed  in  three  stages.  In 
the  first  stage,  tne  dimensions  of  the  piece  appear  to  change  insignificantly,  if 
at  all,  but  tne  mass  drops  as  the  outer  portions  pyrolyze.     In  tne  final  stage,  the 
piece  remaining  is  entirely  reduced  to  char  and  mass  loss  proceeds  from  shrinkage 
of  tne  size  of  the  piece  as  it  bums  away  by  glowing  combustion.     Between  these  two 
stages  a  transition  phase  exists  during  wnich  a  char  layer  is  being  formed  and  being 
eroded  simultaneously.     Comprehensive  models  of  this  process  have  been  attempted 
by  several  workers  (see,  e.g.,  Kanury  1974;  Kung  1971)  but  they  are  so  complicated 
that  tne  level  of  refinement  of  the  rest  of  this  analysis  precludes  their  use  here. 
Hie  simplistic  model  employed  nere  should  suffice  for  order-of-magnitude  computation 

As  shown  on  the  graph  of  figure  C-1,  the  regression  coefficient  value  obtained 

was 

b  =  K  =  0.0064.  (CIO) 


This  number  was  used  in  deriving  the  maximum  range  formulas  graphed  in  the  text. 
Revision  of  this  number  by  including  more  experimental  data  can  easily  be  reflected 
in  tne  results,  as  is  shown  in  appendix  F. 
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Records  had  been  discarded  for  lack  of  storage  space. 
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APPENDIX  D:  LOFTING  OF  FIREBRANDS  BY  FLAME 

AND  BUOYANT  PLUME 


In  tiiis  appendix  a  model  is  developed  for  the  process  by  whicli  potential  fire- 
brands are  lofted  into  tne  ouoyant  plume  aoove  a  steady  flame.     Tne  flame  is,  of 
course,  from  a  burning  tree  or  group  of  trees,  so  is  not  truly  steady,  but  transient. 
Tne  flame/buoyant  plume  structure  is  treated  as  a  steady  flow  field  out  is  permitted 
to  collapse,  from  the  bottom  upward,  when  the  foliage  of  the  burning  tree(s}  is 
consumed.     The  potential  firebrand  is  represented  in  this  model  as  an  insert  cylin- 
der with  specific  gravity  typical  of  weathered  wood. 

Tne  analysis  proceeds  in  three  phases.     First  an  approximate  formula  is  deve- 
loped for  tne  velocity  of  a  firebrand  in  the  vertical  flow  field.     This  formula  is 
then  applied  to  derive  tne  height  of  the  particle  as  a  function  of  time,  assuming 
that  the  flow  field  is  truly  steady.     Tne  demise  of  the  flow  field  is  then  approxi- 
mated by  using  the  steady  state  gas  velocity  versus  height  equations  to  predict  the 
progress  of  the  last  oit  of  flame  product  upward  through  the  flow  structure.  IVith 
the  passage  of  tnis  imaginary  surface,  the  flow  is  presumed  to  be  terminated  as  a 
conerent  structure,  leaving  any  particles  "overtaken"  to  be  transported  by  the  wind 
field  from  tne  height  at  which  they  are  when  overtaken. 

I.     Veloaity  of  an  Inert  Farticle  in  a  Constant  Vertical  Flew  Field 

The  lifting  of  a  potential  firebrand  by  the  flow  field  within  and  above  a 
burning  tree  caii  oe  described  approximately  by  considering  an  idealization  of  the 
process.     The  particle  starts  presumably  from  rest  and  is  accelerated  by  aerody- 
namic drag  upward  against  the  pull  of  gravity.     As  will  be  shown  below,  particles 
that  have  terminal  fall  velocities  not  greatly  different  from  the  fluid  velocity 
are  accelerated  rapidly  (_in  a  short  distance  as  well  as  in  a  short  time)  to  the 
speed  at  which  the  force  of  gravity  is  just  balanced  by  the  drag  induced  by  its 
velocity  relative  to  the  gas.     This  fact  leads  to  the  useful  approximation  that 
the  particle  is  borne  by  the  (variable)  flow  field  at  such  a  speed  as  to  produce 
the  drag/ gravity  balance  everywhere.     So  the  first  step  in  the  model  development 
is  to  describe  the  acceleration  of  an  inert  particle  from  rest  by  a  constant  flow 
field.     The  following  variables  are  used  in  this  description: 


m  = 

mass  of  inert  particle 

V  = 

vertical  velocity  of  particle 

t  = 

time 

P  = 

gas  density 

S  = 

drag  coefficient  of  particle 

A  = 

area  presented  to  flow  by  particle 

w  = 

gas  velocity  (^vertical) 

g  = 

acceleration  of  gravity. 

Tne  equation  of  motion  can  be  written  as 


"  dt  "  2"  ^S^*^^  "  ~  ^^^^ 

and  put  in  dimensionless  form  by  defining  the  following  variables; 
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1 12 

=  terminal  fall  velocity  in  still  fluid  =  (2mg/pCpA)  '  (D2) 
e  =  gt/v^  (D3) 
s*  =  (w  -  v)/v^  =  normalized  "slip  velocity"  (D4) 
ds*  2 

The  slip  velocity  history  given  by  integrating  (DS)  is 

s*  =  1  +  2/((^  _  exp(2e)   -  l)  (D6) 

^  o  / 

when  tlie  particle  starts  from  rest  at  time  t  =  0.  This  equation  shows  that  the 
ultimate  slip  velocity,  v^,  is  approached  closely  within  a  very  brief  period  of 

time.     For  exanple,  if  v^  is  10  ft/sec,  then  in  one  second  6  has  a  value  of  3.2 

and,  even  if  w  is  very  much  greater  than  v^,  the  slip  velocity  has  fallen  to 

within  0.3  percent  of  its  "final"  value.     The  distance  covered  in  that  time  may 
be  large,  however,  so  (D6)  must  be  integrated  again  to  define  the  spatial  history 
of  tne  particle  velocity.     Denoting  particle  height  by  z,  note  that: 

dz 

—  =  v  =  w  -  vs*.  (D7) 
dt  o 

In  terms  of  dimension  less  variables 

z*  =  gz/v^  (D8) 

w*  =  w/v^  (D9) 

and  using  (06)  for  s*  we  have 

^  -  w*  -  1  .  2(w*  -  l)expC-29) 

de  Cw*  +  1)  -  Cw*  -  l)exp(-20)'  ^^^^^ 

The  last  term  in  CDIO)  is  a  perfect  differential,  so 

z*  =  (w*  -  1)6  -  ln(Cw*  +  1  -   (w*  -  l)exp(-2e))/2)  (Dll) 

if  z  is  zero  at  time  t  =  0.     Recasting  (D6)  in  terms  of  the  ratio  of  particle 
velocity  to  ultimate  particle  velocity, 

V*  =  v/(w  -  v^)  (D12) 

we  can  write  the  velocity  and  height  equations  as : 

V*  =  1  -  2/((w*  +  l)expC2e)  -  (w*  -  1))  (D13) 
z*  =   Cw*  +  1)6  +  ln(l  -  V*).  (D14) 
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These  equations  are  graphed  in  figure  D-1  as  a  crossplot  of  z*  as  a  function  of 
V*.     For  typical  firebrand  particles  we  are  interested  in  "terminal  velocity"  values, 
v^,  of  about  10  feet  per  second  (Clements  1977).     Figure  D-2  shows  actual  vertical 

distances  traveled  by  particles  accelerated  from  rest  in  constant  flow  by  the  time 
they  have  achieved  90  percent  of  their  ultimate  velocities.     This  figure  shows  that 
particles  of  the  type  considered  good  candidates  for  firebrands  achieve  a  "steady" 
slip  velocity  that  oalances  the  force  of  gravity  in  a  very  short  distance  in  the 
flow  conditions  typical  of  the  flames  from  burning  trees   (see  appendixes  A,  B) . 
Thus  the  approximation  is  made: 

The  particle  everywhere  travels  at  the  velocity  which  balances  drag  and 
particle  weight. 

Tliis  approximation  will  overestimate  the  firebrand-loft ing  capability  of  any 
given  flow  structure,  but  not  by  a  substantial  amount.     For  particles  traveling  in 
spatially  variable  flows,  this  approximation  greatly  simplifies  the  analysis.  The 
variable  flow  field  of  the  flame/plume  structure  is  addressed  below. 


FRACTION  OF  ULTIMATE  VELOCITY  AHAINED 


Figure  D-1. — Distance  traveled  by  particle  accelerated  upward  from  rest  by  constant- 
velocity  flow. 
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PARTICLE  TERMINAL  VELOCITY 

Figure  D-2. — Distance  vequived  for  paTtiole  to  accelerate  upward  to  90  percent  of 
ultima:te  velocity ^  starting  from  restj  in  constant  fluid  flow. 


II.    Particle  Trajectories  in  Flame /Plume  Flow  Fields 

The  motion  of  an   inert   particle  entrained  in  the  flow  through  the  upper  part 
of  a  steady  flame  is  approximated  here  by  the  equation 

|-=w-v^  (D15) 

where  the  fluid  velocity,  w,  is  a  function  of  height,  z,  but        is  assumed  to  be 

constant.     TJiis  implies  that  the  density  of  the  flame  gas  is  constant.     The  varia- 
tion of  the  density  is  minor  (appendix  A,  eq.  A49)  and  influences  v^  only  as  the 

square  root  of  density,  so  use  of  the  flame-tip  value  is  justified  here.  The  varia- 
tion of  w  with  z  is  adequately  described  by  the  approximation 

w  =  WpCz/Zp)^"^^  (D16) 

where  the  subscript  F  denotes  flame  tip  value. 

Equation   (D15)  is  readily  integrated  to  yield  the  time,  t„,  required  for  the 

r 

particle  to  travel  from  its  initial  height,  z^ ,  to  the  flame  tip: 
^c^n  1  /o  r         1   -   V  /w 

^    =  1  -  (z  /z  ^  _o  J   r   o_F  1 

2z^  o     F  w„       1  ,     ,     .  1/2  ,  I 

F  F       L(z^/Zp)         -  v^/WpJ 
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By  using  the  expression  for  Wp  developed  in  appendix  A  and  by  expressing  in 

terms  of  the  diameter  of  a  cylinder  the  pertinent  times  and  physical  sizes  can 
be  made  more  apparent : 

1/2 

v^Wp  =  B(D/Zp}  (D18} 

where 

B  =  (l^f^y  fl  -  P/P  J  (D19) 


and 


=  drag  coefficient  =  1.2  for  cylinder  in  cross  flow 

=  density  of  cylinder 
p     =  density  of  ambient  air 
Pp  =  density  of  flame  gases  at  flame  tip. 


10 


(D20) 


Tne  parameter  group,  2z  /w      which  divides  the  time  in  equation   (D17)  is  the 

r  r 

time  required  for  fluid  to  travel  from  the  base  to  the  tip  of  the  flame.  This  time 
is  a  characteristic  scale  time  for  the  flow  and  is  used  frequently  in  this  analysis 
to  normalize  travel  time. 

Using  representative  values  for  the  densities  in  expression  CD19) 

p^  =  0.3  gm/cc     (19  Ib/ft"^)  (D21) 
p     =  1.2  X  10"^  gm/cc     (0.075  Ib/ft"^)  (D22) 
Pp/p^  =  3/8.  (D23} 

The  value  of  B  is  found  to  De  about  40.  This  value  was  used  to  construct  figure 
D-3.     This  figure  shows  the  travel  time  from  initial  height   (z^)  to  flame  tip  as 

a  function  of  D/Zp.     The  time  becomes  indefinitely  large  when  the  fluid  velocity 

at  the  initial  height  approaches  the  terminal  velocity  of  the  cylinder  in  the  flame 
gas.     But  for  those  cylinders  that  start  at  half  the  flame  height  or  greater,  the 
travel  time  is  generally  less  tnan  1.5  unless  the  diameter  is  nearly  equal  to  the 
value  that  will  not  permit  the  particle  to  be  lifted  by  the  flow. 

In  the  ''transition  zone"  between  the  tip  of  the  flame  and  the  buoyant  plume, 
the  particle's  terminal  velocity  decreases  because  the  fluid  density  increases, 
but  the  fluid  flow  velocity  is  essentially  constant: 


■'•'^Most    objects,  unless  aerodynamical ly  shaped  for  stable  flight,  will  tend  to 
assume  the  orientation  that  maximizes  drag  when  in  unconstrained  flight.     The  drag 
coefficient  for  cylinders   (or  flat  disks)  to  be  used  here  is  about  1.2. 
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where 

1/2 

and        is  the  particle's  terminal  velocity  in  fluid  with  a  density  equal  to  that 

at  the  flame  tip.     The  travel  time,  t   ,  through  the  transition  zone  (from  z  =  z     to  z  = 

t  r 

1.4z„)  is  given  by  the  integral  of  CD24) .     Expressed  in  dimensionless  form  as  before 

F 

this  time  is : 


0  5  10  15  20  25  30 

lO'X  CYLINDER  DIAMETER  / STEADY  FLAME  HEIGHT 

Figure  D-3. — Approximate  travel  time  for  wood  cylinder  lifted  hy  flame. 
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Since 


V  /w    =  B(D/z  ) 
or  r 


1/2 


(D28) 


and 


(1.4)^  =  2, 


(D29) 


this  can  be  written  as 


,t^/2zp  =  0.2  .  -  B(^)'''ln(l  .  1/ (l 


1/2 


w 


1/2, 


(D30) 


This  expression,  like  equation  (D17) ,  gives  travel  times  near  unity  over  the 
pertinent  range  of  D/Zp.     Figure  D-4  is  a  graph  of  equation   (D30) . 

To  the  travel  times  given  in  figures  D-3  and  D-4  must  be  added  a  travel  time 
within  the  plume.     This  last  time  will  be  a  function  of  the  height  to  which  the 
particle  is  lifted.     Tlien  by  comparing  the  time  it  takes  the  particle  to  reach  a 
particular  height  with  the  time  the  flow  structure  will  persist  to  that  height  we 
can  derive  tne  maximum  achievable  height  as  a  function  of  particle  size. 
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Figure  d-4. — Approximate  travel  time  for  uood  cylinder  through  transition 
zone  hetideen  tip  of  flame  and  buoyant  plum. 
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The  flight  of  the  particle  in  the  buoyant  plume  is  approximated  here  by 

^  =  w  -  V  (D31) 
dt  op  ^ 

where,  from  appendix  B,  the  gas  velocity  is  approximated  by 

w  =  Wp(a/(b  +  z/Zp))^^^  (D32) 

with 

a  =  5.963  (D33) 
b  =  a  -  1.4  =  4.563.  (D34) 

The  particle  terminal  velocity  in  the  plume,  v     ,  is  related  to  that  in  the  flame, 

V  ,  by 
o '  ^ 

1/2 

V      =  O.Sv    =  0.8Bw„(D/zJ^    .  (035) 
op  o  F        F  ^ 


Again  normalizing  the  particle  height  by  the  steady  flame  height  and  the  time  by 
the  characteristic  gas  travel  time  through  the  flame,  the  time,  t^,  required  for  a 

given  particle  to  travel  from  the  base  of  the  plume  (1.4Zp)  to  the  height  z  can  be 
written  as 

.'Z/Z 

Wpt 

^dx.  (D36) 


t        .     f      F/  1/2  \  -, 

f-J  ((b-^T)  -0.8v„/w^)-^ 

•'a  _  n 


2z 

t 

•'a-D 

Tne  result  of  the  indicated  integration  in  (D36)  can  be  expressed  conveniently  in 
terms  of  the  variaole  r,  where 

1  /2 

r  =  CCb  +  z/Zp)/a)       .  (D37) 
In  terms  of  r,  equation  (D36)  gives 

2 

w  t  r/l-0.8vw      \  V.  v\„  "1 

F/v\^  o     F  F  F  ' 

0.8-^ 

Wp 

Again,  using  (D18)  witn  a  value  of  40  for  B,  this  formula  can  be  used  to  calculate 
particle  travel  time  in  the  plume.     By  combining  the  times  given  by  equations  (D38) , 
(D30) ,  and  (D17)  one  obtains  the  total  travel  time  from  initial  particle  height, 
z^,  to  the  "final"  height,  z. 
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Tlie  contrioution  to  this  total  time  from  tne  variation  in  initial  particle 
heignt  (figure  D-S)   complicates  the  calculation  and  presentation  of  the  total  time 
of  travel.     The  refinement  of  a  continuously  variable  initial  height  is  clearly 
not  warranted  ay  tne  precision  of  the  models  presented  here,  either.     So  as  a 
compromise  between  unwarranted  generality  and  oversimplification,  four  representa- 
tive values  of  ^^/-p  \^'ere  diosen  for  display  and  further  computation.     Because  the 

initial  neight  of  a  firebrand  particle  would  be  no  greater  than  the  height  of  a 
tree  contriouting  to  the  flame  and  because  it  has  been  assumed  that  the  base  of 
the  flame  is  at  one-half  the  tree  heiglit,  the  following  representative  values 
were  chosen; 


Tree  Height /Flame  Height   Change)        <0.5        .5-1.0        1.0-1.5  >1.5 

Particle  Initial  Height/Fl ame 

Height   (Range)  <  0.25      .25-. 5  .5-. 75  >.75 

Nominal  Value  ofz/z^  0.2  0.4  0.6  0.8 

o  t 


These  values  probably  typify  the  range  of  conditions  to  be  encountered  when  single 
trees  or  groups  of  trees  "torch  out." 


Using  the  nominal  values  for  z  /z„  in  the  tabulation  above,  figures  D-5 

C  O        F  'to 

through  D-8  were  constructed  by  adding  the  travel  times  from  equations  (D17) , 
(D30] ,     and  (D38) .     Because  the  phase  of  travel  from  initial  height  to  flame 
tip  is  a  strong  function  of  particle  diameter  for  small  values  of  zjz^  (i.e., 

for  tall  flames)  the  "resolution"  in  terms  of  D/Zp  varies  from  one  figure  to 
another. 

These  figures  provide  the  basis  for  a  graphical  solution  of  the  firebrand 
lifting  proolem.     Tne  trajectories  described  by  these  four  figures  require  that 
the  gas  flow  exist  over  the  height  and  time  span  indicated.     So  where  the  line 
describing  the  trajectory  of  the  last  bit  of  flame  gas  crosses  one  of  tlie  par- 
ticle trajectories,  it  marks  the  maximum  height  achievable  by  a  particle  of  the 
size  associated  with  the  intersected  trajectory.     Tl"ie  time  required  for  the  col- 
lapse of  tlie  flame/plume  structure,  in  dimensionless  units,  is  readily  available 
from  the  equations  already  presented  (D17,  D50 ,  and  D38)  .     One  merely  takes  the 
limiting  forms  of  these  equations  when  the  particle  diameter  approaches  zero. 
Doing  so  gives  tne  following  results : 


Time  to  travel  from  flame  base  (z^  =  0)  to  flame  tip 


lim  I 

V  ->  0 
o 


1/2 


In 


1 


V  /w 
o  F 


c     I     ^1/2  / 


t^    (eq.  D17) 


=  1 


(D39) 


Time  to  travel  through  "transition  zone"  =  t^   (eq.  D30) 


0.2 


n    1/2.  ^    1/2      ,  ^  1/2 

B(l-)       (,  .B(5-)      l„(l  .  1/(1  -  B(f-)  )))} 


=  0.2 


(D40) 


53 


Time  to  travel  from  plume  base  to  height  z  =  t^.  For  this  expression  it  is 
simpler  to  approach  the  limit  in  equation  (D36)  before  integrating: 

o  •'a-D  "^a-b 


(D42) 

Tne  total  time  for  which  the  flow  structure  exists  to  height  z,  then  is  given 

3/2 


by  t^,  where 


=  ^o  ^  ^  ^  h  '  ^  =        "  "  Ti  )       -  1^  (D43) 

and  t^  is  tne  period  of  steady  burning  of  the  tree  crown  (s)  normalized  by  the 
characteristic  time  2Zp/Wp.     Equation  (D43)  is  plotted-^  ^  in  figure  D-9  on  the  same 

scale  as  figures  D-5  tnrough  D-8,  using  various  values  for  t^.     By  overlaying 

figures  D-5  tnrougii  D-8  on  this  graph,  one  can  derive  the  maximum  particle  height 

as  a  function  of  particle  size  aiid  steady  flame  duration,  t  . 

o 

Alternatively,  one  may  insert  the  values  of  t    and  z/z    obtained  from  the 

1  r 

particle  trajectory  figures  into  equation  (043)  and  derive  the  required  flame  dura- 
tion, t^.     Using  intermediate  values  not  necessarily  displayed  in  figures  D-5 

tnrough  D-8,  figures  D-10  tnrough  D-13  were  derived  in  this  way.     Tnese  figures 
permit  the  final  step  in  the  long  procedure  outlined  here--fixing  the  maximum 
height  of  viable  firebrands. 

Overlaid  on  tnese  figures  is  the  dashed  line  relating  the  initial  diameter  of 
a  wood  particle  and  tne  maximum  height  through  which  it  may  fall  and  yet  be  burning 
(glowing)  on  impact.     This  relationship,  based  on  the  wind  model  presented  in 
appendix  E  and  trajectory  analysis  described  in  appendix  F,  is: 

max(z)  =  0.39  x  lO^D.  (044) 

The  initial  particle  iieight  must  be  added  to  the  height  increment  gained  in  the 
flame  and  plume,  so  the  dashed  line  must  intercept  the  solid  lines  at  the  points 
where 

z/z„  =  0.39  X  10^Ci^/z„)  -  z  /z„.  (D45) 
r  r  O  r 

The  final  result,  a  plot  of  tne  maximum  viable  particle  height  versus  flame  duration 
for  the  four  initial  heights,  is  given  as  figure  7  in  the  text. 


^^Tne  plot  is  readily  extended  oelow  z/Zp  =  1;  the  functional  form  of  t^^  for 
z  <  z     is  a  square  root. 
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Figure  D-5. — Total  travel  time  for  wood 
cylinder  lifted  by,  flame  into  buoyant 
plume.     Initial  particle  height  = 
0.2  X  steady  flame  height. 
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Figure  D-6. — Total  travel  time  for  wood 
cylinder  lifted  by  flame  into  buoyant 
plume.     Initial  particle  height  = 


0.4  X  steady  flame  height.  PARTICLE  HEIGHT /STEADY  flame  HEIGHT 
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Figure  D-7 . — Total  travel  time  for  wood 
cylinder  lifted  by  flame  into  buoyant 
plume.     Initial  particle  height  = 
0.6  X  steady  flame  height. 
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Figure  D-8. — Total  travel  time  for  wood 
cylinder  lifted  by  flame  into  buoyant 
plume.     Initial  particle  height  = 
0.8  X  steady  flame  height. 


56 


Figure  D-9. — Total  gas  flow  persistence 
time  for  buoyant  plume  above  flame. 
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Figure  D-10. — Flame  duration  required  to  lift  wood  aijlind£.r  to  various  heights^ 
stoj-ting  from  0.2  x  steady  flame  height. 
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Figure  D-11. — Flame  durcction  required  to  lift  wood  cylinder  to  various  heights^ 
starting  from  0.4  x  steady  flam  height. 
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Figuj'e  D-12. — Flame  duration  required  to  lift  wood  chjlinder  to  various  heights^ 
starting  from  0.6  x  steady  flame  height. 
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FigvLPe  D-13. — Flame  duraticn  veqidred  to  lift  wood  aylinder  to  various  heights ^ 
starting  from  0.8  x  steady  flame  height. 
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APPENDIX  E:  A  CRUDE  MODEL  FOR  SURFACE  WIND 

OVER  ROUGH  TERRAIN 


In  order  to  calculate  the  distance  tiiat  a  firebrand  might  be  carried  by  the 
wind,  it  is  necessary  to  descriue  the  wind  field  in  detail.     Such  a  description  is 
clearly  beyond  the  state  of  tne  art  of  meteorology  when  the  problem  is  stated  in 
general  terms.     That  is,  the  details  of  the  wind  velocity  pattern  over  rough  terrain 
for  arbitrary  atmospheric  conditions  are  not  predictable  at  the  present  time.  But 
since  we  are  concerned  only  with  making  rather  crude  "limiting"  calculations,  it 
is  sufficient  for  the  purpose  at  nand  to  describe  the  general  features  of  the  flow 
field  in  some  typical  situations.     Tlie  crude  model  should  reflect  the  influence  of 
tne  terrain  and  tlie  roughness  of  the  surface  cover  on  the  windspeed,  but  it  is  not 
necessary  that  the  model  oe  forecastaDle  or  even  relataole  to  atmospheric  transport 
processes.     It  need  only  be  descriptive,  and  correct  in  the  sense  and  approximate 
magnitude  of  factors  tiiat  influence  the  wind  pattern. 

Two  idealized  situations  are  used  here  to  typify  the  spotting  problem: 

1.  Tl-ie  terrain  can  be  considered  to  be  flat  in  the  direction  of  the  prevailing 

wind. 

2.  The  wind  blows  directly  across  a  "smooth"  ridge/valley  profile  that  repre- 
sents the  major  terrain  relief  in  the  area  of  concern. 

These  two  situations  obviously  do  not  exiiaust  the  possibilities,  but  until 
the  entire  spot  fire  model  is  validated  to  some  degree  in  these  easily  described 
situations  tliere  seems  little  point  to  complicating  tlie  picture  further, 

Tne  first  case  is  intended  to  represent  not  only  genuinely  flat  or  gentle 
terrain,  but  airflow  up  or  down  drainages.     Tlie  second  case  is  intended  to  repre- 
sent a  typical  fire  control  problem--spotting  over  a  ridgetop  or  valley  floor  fuel 
break. 

For  the  "flat  terrain"  case  we  shall  use  only  a  horizontal  wind  component. 
The  windspeed  near  the  surface,  above  the  canopy  of  a  forested  area,  is  well  de- 
scribed by  a  logarithmic  profile  CSutton  1953;  Marunich  1975}  : 

U  =  CbyK)ln(z/z^)  (El) 

where 

U  =  Horizontal  windspeed  Cohort  term  average) 
z  =  height  above  ground 
z^  =  "roughness  length,"  a  linear  scale  associated  with  tlie 
turbulent  transport  of  shear  stress 
=  "friction  velocity,"  a  cnaracteristic  windspeed 
K  =  von  Karman's  constant,  usually  reckoned  at  0.4. 

Data  presented  by  Marunich  (19  75)  and  Tanner  and  Pelton  (1960)  indicate  that 
the  roughness  length  is  related  to  tree  height  for  wind  over  a  uniform-height 
forest  canopy.     These  data  provide  the  following  values: 

z     =  0.1313H  (E2) 

=  °-l^"lz=1.6H  ^E^^ 
where  H  =  tree  iieigiit. 
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Other  surface  covers  give  rise  to  other  (smaller)  friction  lengths  (Sutton  1953). 
Referring  velocities  to  the  value  just  at  treetop  height,  U„,  gives 


This  is  the  form  used  in  tlie  calculation  of  firebrand  transport  presented  in  the 
text . 

For  the  second  case,  the  matter  is  not  so  simple.     In  this  case  the  velocity  of 
the  wind  in  the  vertical  direction  (W)  must  be  taken  into  account,  as  well  as  the 
variation  of  windspeed  with  height  (z)  and  horizontal  distance  (X).     It  has  often 
been  noted  that  windspeed  is  greater  at  ridgetops  than  at  mi ds lope,  and  greater 
there  than  on  tiie  valley  floor  (Sutton  1953).     Tlie  modeling  of  these  variations 
has  met  with  little  success  (Brown  1974)  and  is  the  subject  of  current  research 
mid  development  effort  (Fosberg  and  others  1976;  Ryan  1977). 

Here  we  wish  to  capture  only  the  gross  features  of  such  variations  in  a  simple 
form.     To  do  this  we  posit  an  oversimplified  physical  picture  of  the  flow  field 
over  the  rough  terrain.     Tne  physical  assumptions  of  this  model  are: 

1.  There  exists  an  inviscid  flow  field  trapped  between  an  isobar  at  (constant) 
height  z^  and  a  friction-dominated  surface  layer  of  height  z„r,(^)- 

2.  The  flow  fields  in  botii  the  inviscid  upper  layer  and  the  friction  layer 
are  two-dimensional.     Tnat  is,  velocity  exists  in  the  X  direction  (U)  and  the  z 
direction  (W)  only. 

3.  The  friction  layer  boundary,  z„  (X) ,  is  a  streamline  below  which  the 

r  R 

horizontal  velocity  is  logarithmic  in  z  with  a  constant  roughness  length,  z  . 
No  stability  correction  is  applied. 

4.  Tlie  flow  satisfies  the  incompressible  continuity  equation  everywhere  and 
air  density  variations  are  neglected. 

5.  Horizontal  velocity  is  independent  of  height  in  the  inviscid  layer.  In 
crude  form,  then,  this  model  mimics  some  major  features  of  Fosberg 's  (19  76)  model. 
In  the  model  presented  by  Fosberg  it  is  suggested  that  the  "lid"  height,  z  , 

should  be  taken  to  be  1.5  to  2.0  km  above  the  mean  height  of  the  terrain.     For  the 
purpose  of  presenting  numerical  results  in  this  paper  we  take  z^  to  be  6,000  ft 
(about  1.8  km)  aoove  the  mean  terrain  height. 

The  flow  field  resulting  from  the  above  assumptions  can  be  derived  in  a 
straightforward  way.     First,  in  the  inviscid  upper  layer  we  have: 


U=  U„ln(z/0. 1313H)/ln(l/0. 1313)  =  0 . 49  3U„ln(z/0 . 1313H) . 


(E4) 


constant  =  V 


CE5) 


From  the  continuity  equation 


9u  m_ 

3X  ^  dz 


(E6) 
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tilis  implies  a  vertical  velocity  distribution 


V(X,z)  =  -[V,z/(z^  -  Zp^) 


2  ^"FR 
dX 


f(X). 


(E7) 


Hie  arbitrary  function  fCX)  is  determined  by  matching  the  upper  boundary  condition 

tiiat  z    is  constant,  a  horizontal  streamline: 
(j 


W(X,z)  = 


FR 


7 


dz 


FR 


dX 


=  U 


-G  ~  ^FR_ 


dz 


FR 


dX 


(ES) 


In  the  friction  layer  oetween  the  terrain  surface,  z^(X) ,  and  the  streamline 
Zpj^(X)  ,  we  have 


UJX)  fZ  -  z  (X) 
U  =  — rr—  In' 


K 


o 


(E9) 


Tne  arbitrary  friction  velocity  allows  the  matching  of  horizontal  windspeed  at 
the  friction  layer  upper  boundary,  with  the  result; 


U  = 


-  'fR^^^ 


In 


In 


Zpj^(X)   -  z^(X) 


(£10) 


Tne  vertical  velocity  field  in  the  friction  layer  is  to  be  found  by  integrating 
the  continuity  equation.     The  result  can  be  abbreviated  as: 


„  =  u  — t  -  z  § 


[ 


In  y  dy,  where 


(Ell) 


F(X)  = 


V  z      -  z 

G  ,    /  FR  T 

In 


"g  -  ^FR^^^ 


o 


(E12) 


Equation  (Ell)  satisfies  the  boundary  condition  at  the  terrain  surface  (z  =  +  z^)  ; 
by  imposing  the  boundary  condition  at  the  top  of  the  friction  layer 


FR 


(E13) 


we  obtain  an  equation  describing  the  dependence  of  z  „  ori  z  .     This  equation  is 

r  R  i 


dz 


FR 


dz 


T 


dX  dX 


-  z  (z„ 
o  G 


^FR^ 


f 


m  y  dy  ^((z^ 


.r)-(^)) 


-1 


(E14) 
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The  dependence  of  Zp^  on        can  best  be  deduced  by  relating  two  intermediate 
variables : 

f  =  (Zpj^  -  z^)/z^  CE15) 
g  =  ^^G  -  ^FR^/^o  ^E16) 

g=  -g^  m  y  dy  ^(^^)^  ^^^^^ 

Carrying  out  the  indicated  differentiation  and  rearranging  terms  this  becomes  an 
easily  integrated  form: 


In 


In 


In  y  dy|  -  ln|ln  f j    =  In  g  +  constant  CE19) 


or 


/  In  y  dy  /g  In  f  =  constant  =  C.  (£20) 

This  last  form  is  most  revealing  if  cast  in  terms  of  the  ratio  of  friction  layer 
thickness  to  total  (friction  plus  inviscid)   layer  thickness 

^FR  ~  ^T 

$  =  ^        =  f/Cf  +  g)  (E21) 

T 

$  =  (C  -  1/U  In  f))/(c  +  1  -  (1  +  -)/ln  f).  (E22) 

Because  both  f  and  g  are  large  numbers  relative  to  unity,  equation  (E22)  provides 
the  extremely  useful  result,  that  $  may  be  considered  to  be  constant.     From  equa- 
tion (E20)  it  is  clear  that  C  ^and  hence  $)  is  a  small  number.     In  fact,  Sutton 
(1953)  cites  data  taken  over  flat  terrain  that  indicate  $  =  0.12.     Brown  (1974) 
"splices"  friction  layer  profiles  to  Ekman  layer  profiles  and  asserts  that  this 
layer  depth  ratio  should  be  of  the  order  of  0.1.     In  numerical  results  presented 
herein  we  have  used  $  =  0.12. 

A  further  simplification  obtains  from  this  result  and  equation  (Ell).  Note 
from  (Ell)  that  the  relationship  between  W  and  U  in  the  friction  layer  varies  be- 
tween the  limits 
dz 

^  =  \]  atz  =  z^+z  (E23) 

dX  T  o 

and 

dz  dz 

W  =  U  =  (1  -  <j,)u  ^      at  z  =  Zpj^.  (E24) 
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Because  the  "error"  is  so  slight  and  the  simplification  is  so  great  in  the  particle 
trajectory  equations,  the  furtlier  approximation  is  made  that  the  streamlines  follow 
the  terrain  contour  in  the  friction  layer: 


dz 


IV  =  U 


T 


dX 


CE25) 


From  equations   (E5) ,   (E8) ,  and  (ElO)  it  is  also  clear  that  the  friction  layer 
equations  for  U  and  W  can  be  extended  into  the  in viscid  zone  for  a  distance  equal 
to  a  large  fraction  of  the  friction  layer  height  before  significant  error  accumu- 
lates.    This  approximation  is  also  invoked  in  deriving  the  firebrand  trajectories 
in  appendix  F. 

The  variation  of  horizontal  wind  velocity  with  distance  (UC^))  is  plotted  in 
figure  E-1  for  terrain  described  by  a  sinusoid: 


z^(X)  =  a  sin  (ttX/X^) 


CE26) 


The  total  terrain  relief,  2a,  is  a  parameter  for  these  curves 
are  given  by  equation  (ElO)  with 


The  curves  plotted 


FR 


-  z^)  =  0.12 


(E27) 


and  z    equal  to  6,000  ft  above  the  average  value  of  z„.     This  plot  shows  that  for 

a  ridge-to-valley  elevation  change  of  2,000  ft,  the  velocity  at  a  constant  height 
above  the  terrain  differs  by  aoout  40  percent  between  ridge  and  valley  locations. 
It  is  this  variation  of  horizontal  windspeed  with  horizontal  distance  that  causes 
the  differences  in  potential  spotting  distances  shown  in  the  graphs  in  the  text. 

To  facilitate  transformation  of  windspeed  measurements  or  forecasts  (to  estab- 
lish the  midslope  reference  value  of  treetop  height  wind),  the  variation  of  hori- 
zontal windspeed  with  terrain  elevation  is  graphed  in  the  text. 


Figure  E-l. — Variation  of  horizontal  wind- 
speed  at  constant  height  above  terrain 
for  flow  across  ridge  and  valley.  Model 
results  for  sinusoid  terrain  profile  and 
uniform  forest  cover ,  100  ft  (SO  m)  trees. 
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APPENDIX  F:  FIREBRAND  TRAJECTORIES 


In  this  appendix  we  combine  the  model  for  the  burning  of  wood  cylinders  [appen- 
dix C)  and  the  model  for  the  windspeed  over  rough  terrain  (appendix  E)  to  calculate 
the  maximum  possible  travel  distance  of  such  particles.     This  is  done  in  two  steps. 
First  we  compute  the  maximum  spot  fire  distance  over  flat  terrain,  then  we  calculate 
a  correction  factor  for  rough  terrain. 


Using  the  approximation  that  the  particle  ever>nvhere  has  a  relative  velocity  in 
the  vertical  direction  equal  to  its  instantaneous  terminal  fall  velocity  (appendix 
D)  we  can  write  the  equations  of  motion  for  the  burning  firebrand  as  follows: 

^  =  U(X,z)  (Fl) 

^  =  W(X,z)  -  v^Ct).  (F2) 

Here  we  use  the  symbols 

X  =  horizontal  distance 
z  =  vertical  distance 

t  =  time  (since  start  of  firebrand  fall) 
U  =  horizontal  windspeed 
W  =  vertical  windspeed 
v^  =  particle  terminal  fall  velocity,  a  function  of  t. 

The  fact  that  the  relative  velocity  of  the  particle  with  respect  to  the  wind 
is  everywhere  equal  to  v^  allows  us  to  obtain  the  time  dependence  of  v^  directly. 

From  appendix  C  we  have  the  model 

■^Cp  D)  =  -Kp  V  (F3) 
dt     s  ^  a  o  ^ 

where 

K  =  0.0064 

=  density  of  wood  cylinder 

D  =  mean  diameter  of  wood  cylinder 

p    =  density  of  air 
a 

v^  =  relative  velocity  of  cylinder  and  air  Qiere  equal  to  terminal  velocity). 

But  since  the  terminal  velocity  is  given  by 

=  C^gP3D/2C^p^)^/^  (F4) 

where 


g  =  acceleration  of  gravity 
=  drag  coefficient  (1.2) 

equation  (F3)  can  be  written  in  terms  of  v 
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p  D  =  2C-,p  V  /irg 
s  Da  o  ^ 


dv 


dv 


■Kp  v 
a  o 


This  leads  immediately  to  a  most  useful  result 
v^Ct)  =  v^CO)(l  -  K7Tgt/4CpV^(0))  . 


(F5) 
(F6) 
(F7) 

(F8) 


In  other  words,  the  particle  slows  down  its  relative  fall  velocity  linearly  in  time 
During  the  time  it  bums  out,  it  moves  a  distance  Az  relative  to  the  air  in  which 
it  is  borne  along,  where 


Az  =  ^  V  (0) 
2     o  ^ 


KTTg 


(P,D)^/Kp^ 


(F9) 


and 


(p  D)  =  the  initial  value  of  (p  D) 
^  s    o  s 


Because,  according  to  the  wind  model  derived  in  appendix  E,  streamlines  near 
the  surface  are  parallel  to  the  terrain  surface,  this  "fall  distance"  is  a  constant 
of  the  trajectory  and  is  equal  to  the  initial  height  of  the  particle  if  the  par- 
ticle bums  out  on  retuming  to  the  surface.     This  fact  was  used  in  appendix  D  to 
identify  those  particles  that  have  the  greatest  potential  spot  distance. 


Using  the  shorthand  notation 


T    =  4CpV^C0)/KTTg 


equation  (F2)  can  be  written 

^  =  W(X,z)   -  v^(0)(l  -  t/x) 


(FlO) 


CFll) 


Flat  Tevvain  Spot  Distance 

Equations  (Fl)  and  (Fll) ,  with  U  given  by  a  logarithmic  function  of  z  and  W 
equal  to  zero,  describe  the  flight  of  a  firebrand  over  flat  terrain: 

^=  (U,/K)ln(z/Z^)  (F12) 
^  =  -v^(0}(l  -  t/T).  (F13) 
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The  logarithmic  profile  does  not  accurately  describe  windspeed  variation  below 
treetop  height.     Although  such  windspeed  variations  have  been  measured  many  times 
(Baynton  1965;  Bergen  1971;  Pons  1940;  Gisbome  1941;  Leonard  and  Federer  1973; 
Martin  1971;  Marunich  1975;  Oliver  1971;  Raynor  1971),  no  universal  rule  is  appli- 
cable because  of  the  many  factors  that  can  influence  the  v/ind  field  structure 
within  and  below  the  crown  layer.     One  approximation  for  the  purpose  at  hand 
would  be  to  use  a  constant  windspeed  below  treetop  height.     This  would  give  a 
maximum  spot  fire  distance,  but  perhaps  one  unrealistically  large.     Another  ap- 
proach would  be  to  consider  the  crown  layer  as  a  mechanical  obstacle  to  the  flight 
of  a  firebrand  and  set  the  windspeed  to  zero  below  treetop  height.     Such  an  appro- 
ximation would  tend  to  give  a  lower  bound  to  the  spot  fire  distance.     A  compromise 
that  interpolates  between  these  extremes  is  to  extend  the  logarithmic  profile  below 
treetop  height  all  the  way  to  zero  velocity.     This  is  the  alternative  used  here. 
The  effect  of  this  assumption  is  assessed  later  in  this  appendix. 

Denoting  treetop  height  by  H  and  the  horizontal  windspeed  at  that  height  by 
U^,  equation  (FIZ)   can  be  written  as 

^  =  U^lnCz/z^3/lnCH/z^) .  (F14) 

The  friction  length,  z^,  is  set  to  0.1313H  for  the  calculation  of  numerical  results. 

The  treetop  height,  and  hence  z   ,  is  considered  to  be  constant  over  the  range  of 
particle  travel. 

Equation  (^13)   can  be  integrated  to  give  the  particle  height  as  a  function 
of  time: 

zCt}  =  zCO)   -  v^C03T(t/T  -  jit/xf).  (F153 

Because  we  seek  the  maximum  spot  fire  distance,  the  particle  should  just  be  con- 
sumed at  the  time  of  contact  with  the  ground,  or 

zCx)  =  0  (F16} 
and  so 

zCt)  =  zCO}(l  -  t/T)^  (F17) 

Using  this  result  in  equation  i¥15)   and  dividing  equation  (F14)  by  the  resulting 
expression  gives  the  equation  for  the  particle's  distance  of  travel  (X)  as  a 
function  of  its  height   C^) : 

-(v^){^f''^^-/^o^^^-Wz^).  (F18) 

This  equation  is  readily  integrated  to  give  the  maximum  horizontal  distance  X*.  Ivhen 
normalized  by  the  tree  height  the  result  is: 
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Since 


ln(H/z^)  =  -ln(0.  151-3)  =  2.03,  (F20) 
this  expression  is  very  well  approximated  by  the  form 

If  the  integral  on  z  had  been  terminated  at  H  rather  than  at        (zero  windspeed 

below  treetop  height)  the  term  /O . 1313  would  be  missing  from  this  expression . 
And  if  the  windspeed  were  set  to  a  constant  for  z  less  than  H,  the  term  /0.1313 
would  be  replaced  by  unity.     Figure  F-1  shows  how  the  results  would  change,  using 
the  constant  windspeed  (maximum  range)  assumption  as  a  reference  condition  for 
comparison . 

Using  equation  (F9)   for  v^(O)  in  terms  of  z(0)    (=  Az)   and  using  numerical 
values  for  the  constants  gives  the  equation  plotted  as  a  nomograph  in  the  text: 

1/2  1/2 
X*  =  21.9UJ^)      [0.362  .  [fl)        1  ln(ii°I)}.  ^  (F22) 

This  result  is  more  cogent  than  it  may  appear  at  first,  since  a  multiplying  factor 
can  extend  the  validity  of  the  expression  to  uneven  terrain. 

Uneven  Terrain  Spot  Distance 

Tlie  flight  of  a  firebrand  over  uneven  terrain  is  described  by  equations  (Fl) 
and  (F16)  when  the  wind  field  is  related  to  the  terrain  surface.     Using  the  wind 
model  described  in  appendix  E  and  especially  the  approximation  that  the  streamlines 
parallel  the  terrain,  we  can  write  the  equation  for  particle  height  as  a  function 
of  time : 

f  =   ar  -  ^(t)  (F23] 

where  z^(X)  =  terrain  height  at  location  X. 

Because  the  particle  travels  horizontally  at  speed  U,  the  height  of  the  par- 
ticle above  ground  obeys  the  simple  equation 

^(z  -  z^)  =  -v^(t)  =  -v^(0)(l  -  t/T).  (F24) 

The  simplification  made  possible  by  this  result  is  enormous.     Because  the 
horizontal  windspeed  is  related  to  the  height  above  the  local  terrain:^- 

v„/ln 


=  U(X,z)  =  ~  ^  ^  In^^-y-^]  (F25) 


See  appendix  E  for  definition  of  symbols  9,  z  ,  and  v 
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and  since 

z  -   2^  =   (z(J3    -   2^(X  =  0})(1  -  t/x)^  =   2*^(1   -   t/T)^  (F26) 

equation  (F23)  is  essentially  the  same  as  equation  (F14) .     If  the  treetop  height 
windspeed,  U^,  is  considered  to  be  a  function  of  liori2ontal  position,  the  equations 

are  identical.     The  integral  forra  for  tne  maximum  spot  distance  is  obtained  as 
before.     Using  equation  (^26]  in  CF24}   and  dividing  CF25)  by  the  result  yields: 

dX   (v,/v„(0])l.HH/z^).(z°/(z  -  z^))l/-ln((z  -  .^)/zJ 

d(z  -  z^J  (1  -  ^)(Zg  -  z^)ln(Kz^  -  z^J/z^)ln(.H/ZgJ        '  ^>-^'> 

The  integral  of  the  function  of  (2  -  2  )  on  the  right  gives  a  function  of  initial 

0 

particle  height  aDove  tne  ground  [z  J   and  tree  height  tnat  is  identical  to  the 
previous  case.     Hence  we  iiave: 


v^m  —  J  /  0^^  ^  ^ 


.^In  —  J,, 
G      2  0 
o 


CF28) 


This  formula  should,  of  course,  reduce  to  the  flat  terrain  form  simply  by 
setting  2^  equal  to  a  constant.     By  choosing  the  constant  value  to  be  the  spatial 
average  terrain  height,  denoted  by  Y^,  we  derive  the  ratio  formula  for  the  maximum 
spotting  distance  in  the  form: 


X*  ln(?  —  

z     -  2         r  2 

°         dX  =  X*Cflat  terrain).  (F29) 


"  ^T  .    /    "G  ~T 


0  In  c 


o 

The  right  hand  side  of  equation  CF29)  is  the  maximum  spot  fire  distance  over  flat 
terrain  using  the  windB^eed  (as  a  fimation  of  'aeight)  tnat  obtains  at  the  average 
tevToln  height. 

The  integral  in  this  equation  can  be  simplified  by  considering  only  cases  for 

whi  ch 

-  Y^|/Cz^  -  z-^)   «  1;  (F303 

for  these  cases,  since 

(z^  -  z^)/Cz^  -  z"^)  =  1  +  (Zj.  -  z^}/(2^  -        ,  (F51) 

^         o       /  \  o     ^        G  T 
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and  so 


Z     -   Z   \\  z      -  z 

1  +  l/ln((})  — -  jj  J  _  --  dX  =  X*(flat  terrain).  (F33) 


The  variation  of  terrain  height  relative  to  the  average  height  can  be  repre- 
sented with  good  generality  by  a  sinusoid, 

z^  -  Ij,  =  a  sin  Cm(X  +  Xp )  (F34) 


wnere 


2a  =  total  elevation  difference,  ridge  to  valley 

m  =  IT  V  horizontal  distance  from  ridge  to  valley 

X^  =  horizontal  location  of  firebrand  source 

(Xj^  =  0  is  midslope  on  windward  side  of  ridge) 


This  form  leads  to  a  transcendental  equation  for  the  maximum  horizontal  distance: 

X*  +     (1  |cos  (m(X*  +  X  ))  -  cos  (mX  )1  =  X*(flat  terrain)  (F35) 

m [  Zg      z^j  L  i  i  J 

where 

a  =  l/ln(Kz^  -  r^O/z^)  .  (F36) 

Equation  (F^S)  was  solved  for  the  ratio  of  maximum  distances  and  plotted  in  the 
text.     The  four  cases  presented  in  the  text  correspond  to  different  values  of  mX^ . 
The  values  used  are  as  follows; 

 1_  Location  of  firebrand  source 

IT  midslope  on  lee  side  of  ridge 

-tt/2  valley  bottom 

0  midslope  on  windward  side  of  ridge 

IT  ridgetop 

For  the  numerical  results  plotted  in  the  text,  the  value  of  a  was  set  to  0.2 
and  z    -  Y    to  6,000  ft.     These  values,  along  with  cj)  =  0.12  (see  appendix  E)  , 

imply  a  roughness  scale  of  z^  =  4.9  ft.     The  weak  dependence  of  the  resulting 

ratio  of  maximum  ranges  on  the  parameter  a  implies  that  the  numerical  value  of  a 
is  not  critical  in  this  range  of  values.     The  location  of  the  firebrand  source 
is  very  important,  however. 
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Figure  F-l.— Effect  of  different 
windspeed  profiles  on  pre- 
dicted maximum  spot  fire 
distance. 
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Burning  Rate  l4odel  Influence 

Note  that  if  the  firebrand  burning  rate  were  different  from  the  model  used 
here  (i.e.,  if  a  different  value  of  the  burning  rate  regression  parameter  b  of 
appendix  C  were  used)  only  tne  value  of  t  would  be  changed.     The  scaling  relation- 
ships would  be  preserved  between  flat  terrain  and  uneven  terrain  results.  The 
impact  of  a  change  in  the  burning  rate  would  be  to  change  the  numerical  coeffi- 
cient K  in  equation  (F9)   and  hence  the  value  of  the  coefficient  in  equation  (F22) . 
The  number  21.9  in  the  latter  equation  is  inversely  proportional  to  the  square 
root  of  the  burning  rate  coefficient  K,  so  results  for  spot  fire  distance  should 
be  altered  accordingly  to  reflect  a  different  K  value. 
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